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We consider two interacting systems when one is treated classically while the 
other system remains quantum. Consistent dynamics of this coupling has been 
shown to exist, and explored in the context of treating space-time classically. 
Here, we prove that any such hybrid dynamics necessarily results in deco- 
herence of the quantum system, and a breakdown in predictability in the 
classical phase space. We further prove that a trade-off between the rate of this 
decoherence and the degree of diffusion induced in the classical system is a 
general feature of all classical quantum dynamics; long coherence times 
require strong diffusion in phase-space relative to the strength of the coupling. 
Applying the trade-off relation to gravity, we find a relationship between the 
strength of gravitationally-induced decoherence versus diffusion of the metric 
and its conjugate momenta. This provides an experimental signature of the- 
ories in which gravity is fundamentally classical. Bounds on decoherence rates 
arising from current interferometry experiments, combined with precision 
measurements of mass, place significant restrictions on theories where Ein- 
stein’s classical theory of gravity interacts with quantum matter. We find that 
part of the parameter space of such theories are already squeezed out, and 
provide figures of merit which can be used in future mass measurements and 
interference experiments. 


When considering the dynamics of composite quantum systems, there 
are many regimes where one system can be taken to be classical and 
the other quantum-mechanical. For example, in quantum thermo- 
dynamics, we often have a quantum system interacting with a large 
thermal reservoir that can be treated classically, whilst in atomic 
physics it is common to consider the behaviour of quantum atoms in 
the presence of classical electromagnetic fields. Things become more 
complicated when one considers classical-quantum (CQ) dynamics 
where the quantum system back-reacts on the classical system. This is 
particularly relevant in gravity because we would like to study the 
back-reaction of thermal radiation being emitted from black holes in 
space-time, and while the matter fields can be described by quantum 


field theory, we only know how to treat space-time classically. Likewise 
in cosmology, vacuum fluctuations are a quantum effect that we 
believe seeds galaxy formation, while the expanding space-time 
they live on can only be treated classically. In addition to the need for 
an effective theory that treats space-time in the classical limit, there 
has long been a debate about whether one should quantise gravity’ ’. 

The prevailing view has been that a quantum-classical coupling is 
inconsistent. Many proposals for such dynamics” are not completely 
positive (CP), meaning they are at best an approximation and fail 
outside a regime of validity’. A map A is completely positive (CP), iff 
1 @ Ais positive. This is the required condition used to derive the GKSL 
equation. If it is violated, the dynamics acting on half of an entangled 
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state, give negative probabilities. The semi-classical Einstein’s 
equation” ’, which replaces the quantum operator corresponding to 
the stress-energy tensor by its expectation value, is another attempt to 
treat the classical limit from an effective point of view, but it is non- 
linear in the state, leading to pathological behaviour if quantum fluc- 
tuations are of comparable magnitude to the stress-energy tensor”. 
This is often the precise regime we would like to understand. 

However, examples of classical-quantum dynamics such as those 
first introduced in refs. 21,22 and studied in refs. 11,23-27 do not suffer 
from such problems and are consistent. More generally, the master 
equation shown in Eq. (4)", is linear in the state space, preserves the 
division of classical degrees of freedom and quantum ones, and is 
completely positive (CP), and preserves normalisation. This ensures 
that probabilities of measurement outcomes remain positive and 
always add to 1. The dynamics is related to the GKSL or Lindblad 
equation’®”’, which for bounded generators of the dynamics, is the 
most general Markovian dynamics for an open quantum system. 
More precisely, we consider dynamics which is autonomous, meaning 
the couplings in the theory do not depend on time. Likewise, Eq. (4) is 
the most general Markovian classical-quantum dynamics with boun- 
ded generators”*°. Sub-classes of this master equation, along with 
measurement and feedback approaches, have been discussed in the 
context of Newtonian models of gravity**”>*!*> and further developed 
into a spatially covariant framework so that Einstein gravity in the ADM 
formalism®’ emerges as a limiting case”. Dynamics which is mani- 
festly diffeomorphism invariant has also been introduced using path 
integral methods***?. 

In this work, we move away from specific realisations of CQ 
dynamics, in order to discuss their common features and the experi- 
mental signatures that follow from this. An early precursor to the 
discussion here is the insight of Didsi”, who added classical noise and 
quantum decoherence to the master equation of ref. 13, and found 
the noise and decoherence trade-off required for the dynamics to 
become completely positive. Here we prove that the phenomena 
found in refs. 11,22,26 are generic features of all CQ dynamics; 
the classical-quantum interaction necessarily induces decoherence on 
the quantum system, and there is a generic trade-off between the 
rate of decoherence and the amount of diffusion in the classical phase 
space. The stronger the interaction between the quantum system and 
the classical one, the greater the trade-off. One cannot have quantum 
systems with long-coherence times without inducing a lot of diffusion 
in the classical system. One can also generalise this result to a trade-off 
between the rate of diffusion and the strength of more general 
couplings to Lindblad operators, with decoherence being a 
special case. 


Results 

Our main result is expressed as Eqs. (26) and (23), which bounds 
the product of diffusion coefficients and Lindblad coupling constants 
in terms of the strength of the CQ-interaction. It is precisely this 
trade-off which allows the theories considered here, to evade the 
no-go arguments of Feynmann’’”, Aharonov’, Eppley and Hannah* 
and _ others'7°51640-48" The essence of arguments against 
quantum-classical interactions is that they would prohibit super- 
positions of quantum systems that source a classical field. Since dif- 
ferent classical fields are perfectly distinguishable in principle, if 
the classical field is in a distinct state for each quantum state in the 
superposition, the classical field could always be used to determine the 
state of the quantum system, causing it to decohere instantly. By 
satisfying the trade-off, the quantum system preserves coherence 
because diffusion of the classical degrees of freedom means that the 
state of the classical field does not determine the state of the quantum 
system”>, Equation (26) and other variants we derive, quantify the 
amount of diffusion required to preserve any amount of coherence. If 
space-time curvature is treated classically, then complete positivity of 


the dynamics means its interaction with quantum fields necessarily 
results in unpredictability and gravitationally induced decoherence. 
This trade-off between the decoherence rate and diffusion pro- 
vides an experimental signature, not only of models of hybrid New- 
tonian dynamics such as refs. 24,33 or post-quantum theories of 
general relativity such as refs. 11,39 but of any theory which treats 
gravity as being fundamentally classical. The metric and their con- 
jugate momenta necessarily diffuse away from what Einstein’s general 
relativity predicts. This experimental signature squeezes 
classical-quantum theories of gravity from both sides: if one has 
shorter decoherence times for superpositions of different mass dis- 
tributions, one necessarily has more diffusion of the metric and con- 
jugate momenta. In the “Methods” subsection “Detecting gravitational 
diffusion” we show that the latter effect causes imprecision in mea- 
surements of mass such as those undertaken in the Cavendish 
experiment*”*' or in measurements of Newton’s constant “Big G"?**. 
The precision at which a mass can be measured in a short time, thus 
provides an upper bound on the amount of gravitational diffusion, as 
quantified by Eq. (42). In the other direction, decoherence experiments 
place a lower bound on the diffusion. Our estimates suggest that 
experimental lower bounds on the coherence time of large 
molecules” °°, combined with gravitational experiments measuring 
the acceleration of small masses*'~’, already place strong restrictions 
on theories where space-time is not quantised. In the section 
“Physical constraints on the classicality of gravity” we show that several 
realisations of CQ-gravity are already ruled out, while other realisations 
produce enough diffusion away from general relativity to be detect- 
able by future table-top experiments. Although the absence of such 
deviations from general relativity would not be as direct a confirmation 
of the quantum nature of gravity, such as experiments proposed in 
refs. 64-74 to exhibit entanglement or coherence generated by grav- 
itons, it would effectively rule out any sensible theory that treats 
space-time classically. While confirmation of gravitational diffusion 
would suggest that space-time is fundamentally classical. Experiments 
to detect or bind gravitational diffusion also provide immediate-term 
prospects for probing the quantumness of gravity, while 
entanglement-based experiments will only be feasible in the long term. 
The outline of the remainder of this paper is as follows. In the 
subsection “Classical-quantum dynamics” we review the general form 
of the CQ master equation of classical-quantum systems as derived in 
refs. 11,30. The CQ-map can be represented in a manner akin to the 
Kraus representation’® for quantum maps, with conditions for it to be 
completely positive and trace preserving (CPTP). We can perform a 
short time moment expansion of the CQ-map taking states at some 
initial time, to states at a later time. This gives us the CQ version of the 
Kramers—Moyal expansion’®*”’, presented in the subsection “The CQ 
Kramers-Moyal expansion”. The physical meaning of the moments is 
given in the subsection “Physical interpretation of the moments”. Our 
main result is presented in the section “A trade-off between deco- 
herence and diffusion”, where we show that there is a general trade-off 
between decoherence of the quantum system and diffusion in the 
classical system. We generalise the trade-off to the case of fields in the 
subsection “Trade-off in the presence of fields” and in the subsection 
“Physical constraints on the classicality of gravity”, we apply the 
inequality in the gravitational setting. The positivity constraints mean 
that the considerations do not depend on the specifics of the theory, 
only that it treats gravity classically, and is time-local. This allows us to 
discuss some of the observational implications of this result and we 
comment on the relevant figures of merit required in interference and 
precision mass measurements in order to constrain theories of gravity, 
as they are not always readily available in published reports. In addition 
to table-top constraints, we consider those due to cosmological 
observations. We then conclude with a discussion of our results. The 
“Methods” section collects or previews a number of technical results. 
Since this paper appeared on the arXiv, we have found that when the 


Nature Communications | (2023)14:7910 


Article 


https://doi.org/10.1038/s41467-023-43348-2 


decoherence-diffusion trade-off is saturated, there are two important 
consequences. The first is that in the continuous class of dynamics, the 
quantum state remains pure, conditioned on the classical trajectory’®. 
The second is that in the path integral formulation, one can show that 
the dynamics are completely positive from the path integral alone”. 
For a generic path integral of Feynman-Vernon form”’, one typically 
only knows that the dynamics are completely positive if it was derived 
from a CPTP master equation. 


Classical-quantum dynamics 

Let us first review the general map and master equation governing 
classical-quantum dynamics. The classical degrees of freedom are 
described by a differential manifold 4 and we shall generically denote 
elements of the classical space by z. For example, we could take the 
classical degrees of freedom to be position and momenta in which case 
M-=R? and z= (q,p). The quantum degrees of freedom are described 
by a Hilbert space H. Given the Hilbert space, we denote the set of 
positive semi-definite operators with trace at most unity as S.,(H). 
Then the CQ object defining the state of the CQ system at a given time 
is a map 9: M-— S.,(H) subject to a normalisation constraint 
J492zT[e]=1. To put it differently, we associate to each classical 
degree of freedom an un-normalised density operator, @(z), such that 
Tr,,[o] = p(z) 2 0 is a normalised probability distribution over the clas- 
sical degrees of freedom and f,,dzQ(z) is a normalised density 
operator on H. An example of such a CQ-state is the CQ qubit depicted 
as a 2 x 2 matrix over phase space”®. More generally, we can define any 
CQ operator f(z) which lives in the fibre bundle with base space M and 
fibre H. 

Just as the Lindblad equation is the most general evolution law 
that maps density matrices to density matrices, we can ask what is the 
most general evolution law that preserves the quantum-classical state- 
space. Any such dynamics, if it is to preserve probabilities, must be 
completely positive, norm-preserving, and linear in the CQ-state. That 
dynamics must be linear can be seen as follows: if someone prepares a 
system in one of two states do or 0, depending on the value of a coin 
toss (|0)(O| with probability p,|1)(1| with probability 1-p), then the 
evolution £ of the system must satisfy p|0)(0| @ Log + (1 — p)|1) (1) @ 
Lo,=L(p09+(1—p)o,) otherwise the system evolves differently 
depending on whether we are aware of the value of the coin toss. A 
violation of linearity further implies that when the system is in state 09 
it evolves differently depending on what state the system would have 
been prepared in, had the coin been |1)(1| instead of |0)(O|. This 
motivates our restriction to linear theories. We will also require the 
map to be Markovian on the combined classical-quantum system, 
which is equivalent to requiring that there is no hidden system that 
acts as a memory. This is natural if the interaction is taken to be fun- 
damental, but is the assumption that one might want to remove if one 
thinks of the hybrid theory as an effective description. We thus take 
these as the minimal requirements that any fundamental 
classical-quantum theory must satisfy if it is to be consistent. 

The most general CQ-dynamics, which maps CQ states onto 
themselves can be written in the form” 


a(z,t+60)= f de’ SN ei2’, 6OL,al2', OL) @) 


yw 


where the L,, is an orthogonal basis of operators and A“”(z|z’, 6t) is 
positive semi-definite for each z,z’. Henceforth, we will adopt 
the Einstein summation convention so that we can drop ¥,, with the 
understanding that equal upper and lower indices are presumed to be 
summed over. The normalisation of probabilities requires 


i, dz (z\z', 6t)LTL, = I. (2) 


The choice of basis L,, is arbitrary, although there may be one 
which allows for unique trajectories”°. Equation (1) can be viewed as a 
generalisation of the Kraus decomposition theorem. 

In the case where the classical degrees of freedom are taken to be 
discrete, Poulin® used the diagonal form of this map to derive the most 
general form of Markovian master equation for bounded operators, 
which is the one introduced in ref. 21. When the classical degrees of 
freedom are taken to live in a continuous configuration space, we need 
to be a little more careful, since e(z) may only be defined in a dis- 
tributional sense; for example, Q(z) =6(z,Z)o(Z). In this case (1) is 
completely positive if the eigenvalues of A“’(z|z’, 6t), A“ (z|z’, 6t), are 
positive so that [ dzdz'P,(z,z')A"(z|z’, 6t) >0O for any vector with 
positive components P,,(z, z’y°. One can derive the CQ master equa- 
tion by performing a short time expansion of Eq. (1) in the case when 
the L,, is bounded”. To do so, we first introduce an arbitrary basis of 
traceless Lindblad operators on the Hilbert space, L,, = {/, L,3. Now, at 
6t=0, we know Eq. (1) is the identity map, which tells us that 
A°°(z|z’, 6t = 0) = 6(z,z’). Looking at the short-time expansion coeffi- 
cients, by Taylor expanding in 6t <1, we can write 


M"(z|z', 6t) = 65660(z, 2) + W#"(z\z')6t + O(6t?). (3) 


By substituting the short-time expansion coefficients into Eq. (1) 
and taking the limit 6t > O we can write the master equation in the form 


OAs 1 sypiteap cane Lapis 
Fes fda! WH ala Wyalz'WL) — SW EILILy a), (4) 


where {, }, is the anti-commutator, and preservation of normalisation 
under the trace and f dz defines 


W#"(z)= / dz’W?"(z'|z). (5) 


We see the CQ master equation is a natural generalisation of the 
Lindblad equation and classical rate equation in the case of 
classical-quantum coupling. We give a more precise interpretation of 
the different terms arising when we perform the Kramers—Moyal 
expansion of the master equation at the end of the section. The posi- 
tivity conditions from Eq. (1) transfer to positivity conditions on the 
master equation via (3). We can write the positivity conditions in an 
illuminating form by writing the short time expansion of the transition 
amplitude A’”(z|z’, 6t), as defined by Eq. (3), in block form 


6(z,z')+ 6tW(z|z") 6twF(zI\z’) 


AY’ (z|z', 6t) = 
lest) btw (z|z’) tw (z\z') 


+O(6t?) (6) 


and the dynamics will be positive if and only if A”’(z|z’, 6t) is a positive 
matrix. It is possible to introduce an arbitrary set of Lindblad operators 
L,, and appropriately redefine the couplings W*”(z|z’) in Eq. (4)". For 
most purposes, we shall work with a set of Lindblad operators that 
includes the identity L,, = (/, La); this is sufficient since any CQ master 
equation is completely positive if and only if it can be brought to the 
form in Eq. (4), where the matrix (6) is positive. 


The CQ Kramers-Moyal expansion 

In order to study the positivity conditions it is first useful to perform a 
moment expansion of the dynamics in a classical-quantum version of 
the Kramers—Moyal expansion as in ref. 11. In classical Markovian 
dynamics, the Kramers—-Moyal expansion relates the master equation 
to the moments of the probability transition amplitude and proves to 
be useful for a multitude of reasons. Firstly, the moments are related to 
observable quantities; for example, the first and second moments of 
the probability transition amplitude characterise the amount of drift 
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and diffusion in the system. This is reviewed in the subsection “Physical 
interpretation of the moments”. Secondly, the positivity conditions on 
the master equation transfer naturally to positivity conditions on the 
moments, which we can then relate to observable quantities. In the 
classical-quantum case, we shall perform a short time moment 
expansion of the transition amplitude A“” (z|z’, 6t) and then show that 
the master equation can be written in terms of these moments. We 
then relate the moments to observational quantities, such as the 
decoherence of the quantum system and the diffusion in the classical 
system. 

We work with the form of the dynamics in Eq. (4), using an arbi- 
trary orthogonal basis of Lindblad operators L,, = {IL}. We take the 
classical degrees of freedom M to be d dimensional, z= (z,, ...zg), and 
we label the components as Z;, i € {1, ...d}. We begin by introducing the 
moments of the transition amplitude W*’(z|z’) appearing in the CQ 
master Eq. (3) 


1 
Dit, @r=7, | WH eRVe-2),..@-Z, — @ 


The subscripts i; € {1, ...d} label the different components of the 
vectors (z — z’). For example, in the case where d= 2 and the classical 
degrees of freedom are position and momenta of a particle, 
Z=(%, 22) =(g,p), then we have 
(Z — 2')=(Z — 24,22 — Z4)=(q — 7’, p — p’). The components are then 
given by (z—2z'),=(q—q) and (z—z’),=(p—p’). Mi, ; (2, 6t) is 
seen to be an nth rank tensor with d” components. 

In terms of the components Dri, oh the short time expansion of 
the transition amplitude A””’(z|z’) is given by 


n 
M” (z|z', 6t) =65.696(z, 2’) + 6t = Dns. i,@) Ca 6(z, 2’) + O(6t?), 
OZ), 


n=0 


(8) 


and the master equation takes the form” 


daz, t) _ n a" 00 
a OO" aa, } (O- 


‘ a . ly é 
— TH(2), a2) ihe. —5DP ZL pla), 9) 


+E See 


pv400 n 


()a(z,0)) 


a} (DY, 1, @Lz. out), 


where we define the Hermitian operator H(z)= fe L, 7 DeFE) 
(which is Hermitian since De = Do OH") We see the first line of Eq. (9) 
describes purely classical dynamics, and is fully described by the 
moments of the identity component of the dynamics A°°(z|z’). The 
second line describes pure quantum Lindbladian evolution described 
by the zeroth moments of the components A®°(z|z’), A"4(z|z’); speci- 
fically the (block) off diagonals, Dp*(2), describe the pure Hamiltonian 
evolution, whilst the components pb (z) describe the dissipative part 
of the pure quantum evolution. Note that the Hamiltonian and Lind- 
blad couplings can depend on the classical degrees of freedom so the 
second line describes the action of the classical system on the quan- 
tum one. The third line contains the non-trivial classical-quantum back- 
reaction, where changes in the distribution over phase space are 
induced and can be accompanied by changes in the quantum state. 


Physical interpretation of the moments 

Let us now briefly review the physical interpretation of the moments 
that will appear in our trade-off relation. In particular, the zeroth 
moment determines the rate of decoherence (and Lindbladian cou- 
pling more generally), the first moment gives the force exerted by the 
quantum system on the classical system, and the second moment 


determines the diffusion of the classical degrees of freedom. For this 
discussion, we shall take the classical degrees of freedom to live in a 
phase space I =(M, @), where @ is the symplectic form. 

Consider the expectation value of any CQ _ operator 
O(z), (O(Z)) := f dzTr[O(z)e] which does not have an explicit time 
dependence. Its evolution law can be determined via Eq. (9) 


OO) 2 ij dztr ov 3 = dzTre 
1 


[0c H(z)] + DEP (ZL ZOR)La — 70 thal, O()} , 
Vv a" t 
Sh) ( gg nag) (ebsO12 8). 


where we have used cyclicity of trace and integration by parts, to bring 
the equation of motion into a form that would enable us to write a CQ 
version of the Heisenberg representation" for a CQ operator. If we are 
interested in the expectation value of phase space variables O(z) =z; 
then Eq. (10) gives 


(10) 


ee f dzDi" Tr}, Lt [LiL yolz, 0) (11) 


with all higher order terms vanishing, and we see that )>,,,,.00D4;( ) 
governs the average rate at which the quantum system moves oe 
classical system through phase space, and with the back-reaction is 
quantified by the Hermitian matrix D{* := (per) . The force of this 
back-reaction is especially apparent if the equations of motion are 
Hamiltonian in the classical limit as in ref. 11. Le. if we define 
H,(2) := AML, and take Dit = wid; n° with w the symplectic form 
and d; the exterior derivative. “Then Eq. (11) is analogous to Hamilton’s 
equations, and the CQ evolution equation after tracing out the 
quantum system has the form of a Liouville’s equation to first order 
and in the classical limit, 


Op(Z, t) 


- (12) 


=(H,, p(Z, 0} + tr ({H)(Z), A(Z)}) + 


with p(Z) := Tr[o(z)]. 
The significance of the second moment is also seen via Eq. (10) to 
be related to the variance of phase space variables 
= (Z;,Z;,) — (2;,) (Zi) 


2,2 rage HEN’ 


OO, 5 pyah aB yt 
“a = 2(D9 Lila) + (22098, Lplia) — (2 


Oh 2, Upla) (13) 


+ (ZDi5, La) - (2, (DEE, LiL). 


In the case when D, z is uncorrelated with z;, and D, ,, uncor- 
related with z, , then the growth of the variance only depends on the 
diffusion coefficient. 

The zeroth moment pi is just the pure Lindbladian couplings. 
The simplest example is the case of a pure decoherence process with a 
single Hermitian Lindblad operator L and decoherence coupling Do. 
Then we can define a basis {|a)} via the eigenvectors of L and 


(a 
and we see that the matrix elements of @ which quantify coherence 
between the states |a), |b) decay exponentially fast with a decay rate of 
Do(L(a)-L(b))?. For a damping/pumping process of a quantum har- 


monic oscillator with Hamiltonian H = qwa'a, Ly =a,L,=a',a the crea- 
tion operator, and D}',D%' the non-zero couplings, then standard 


ee b) = — i(a\(H@),0)16) — 5Do(L(a) — Lib) alelb), 4) 
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calculations***° show that an initial superposition 7; 5|n+m) with n,m 


large and n>m will initially decohere at a rate oF approximately 
(D}* +Ds!)\(m+n)/2, and the state will eventually thermalise to a 
temperature of w/log(D}'/D}'). So in this case, the Lindblad 
couplings not only determine the rate of decoherence but also the 
rate at which energy is pumped into the harmonic oscillator. In the 
next section, we will derive the trade-off between Lindblad couplings 
and the diffusion coefficients. Although we will sometimes refer to this 
as a trade-off between decoherence and diffusion, this terminology 
is only strictly appropriate for pure decoherence processes, while 
more generally, it is a trade-off between Lindblad couplings and 
diffusion coefficients. 


A trade-off between decoherence and diffusion 

In this section, we present our main result by using positivity condi- 
tions to prove the trade-off between decoherence and diffusion seen in 
models such as those of refs. 11,22,26 are in fact a general feature of all 
classical-quantum interactions. We shall also generalise this, and 
derive a trade-off between diffusion and arbitrary Lindbladian cou- 
pling strengths. The trade-off is in relation to the strength of the 
dynamics and is captured by Eqs. (20), (23) and (26). In the subsection 
“Trade off in the presence of fields” we extend the trade-off to the case 
where the classical and quantum degrees of freedom can be fields and 
use this to show that treating the metric as being classical necessarily 
results in diffusion of the gravitational field. 

There are two separate possible sources for the force (or drift) of 
the back-reaction of the quantum system on phase space—it can be 
sourced by either the pee ; components or the Lindbladian components 
De. We shall deal with ‘both sources simultaneously by considering a 
co Cauchy-Schwartz inequality which arises from the positivity of 


Tr / dzdz'N"”(z\z')0,(z,z')p(2’)O'(z,2')| 2 0 (15) 


for any vector of CQ operators O,,. One can verify that this must be 
positive directly from the positivity conditions on A”’(z|z’) and we go 
through the details in the Appendix section “Positivity conditions and 
the trade-off between decoherence and diffusion’. A common choice 
for O, would be the set of operators L,, = {1, L,} appearing in the master 
equation. 

The inequality in Eq. (15) turns out to be especially useful since it 
can be used to define a (pseudo) inner product on a vector of opera- 
tors with components O,, via 


(O,,0,) = / dzdz'Tr [ann Z1z')0 oz’ 03, | (16) 


where ||O|| = \/ (O, O) = 0 due to Eq. (15). Technically this is not positive 
definite, but this shall not be important for our purpose. Taking the 
combination O,,=||O3||?O,, — (O;,02)O, for vectors Oy, O2y, Posi- 
tivity of the norm gives 


|O1|? = |I(!O||7O, — (O,, 02) Oy) || 


y ae oS (17) 
= ||O5||? ({|O41|7||021/? — |(Oy, O2)|7) 20, 
and as long as ||/O,|| #0 we have a Cauchy-Schwartz inequality 
|10,|I7||O||? — (01, 02) |? 2 0. (18) 


We can use Eq. (18) to get a trade-off between the observed dif- 
fusion and decoherence by picking O,,,=6;,L, and O,, =hie= Zily, 
where L,, = {I,L4} are the Lindblad operators appearing in the master 
equation and b'‘ are the components of an arbitrary vector. In this case 
|O,||= f dzTr [Dele ol | and one can verify using CQ Pawula 
theorem®’ that in order to have non-trivial back-reaction on the 


quantum system complete positivity demands that ||O,|| > O, meaning 
the Cauchy-Schwartz inequality in Equation (18) must hold. To reach 
this conclusion one can insert the CQ state into the CQ 
Cauchy-Schwartz inequality and repeat the proof of the Pawula 
theorem*°, which must now hold once averaged over the state. By 
using the short-time moment expansion of A’’(z|z’) defined in Eq. (8) 
and using integration by parts, we then arrive at the observational 
trade-off between decoherence and diffusion 


2 
| dzTr|2" DS ,b/L,(2)L!] / dzTr[DiLqe2)L] 2 


[oz] [ini LZ i] 


, 


(19) 


which must hold for any positive CQ state o(z). Stripping out the b' 
vectors, (19) is equivalent to the matrix positivity condition 


0 < 2(D,)(Dp) — (D?")(D°")", Volz), (20) 


where we define 


(Do) = / dzTr | DL gatz)kp] (DP); = / dzTr [DIL ,o(2)L4), 
(21) 
(D3); - | dzTr/. Le L,o@)L}]. 


Since Eq. (20) holds for all states, the tightest bound is provided 
by the infimum over all states 


0 = inf(2(D2)(Do) — (D’")(D?)"}. (22) 


The quantities (D2) and (Do) appearing in Eq. (20) are related to 
observational quantities. In particular (D2) is the expectation value of 
the amount of classical diffusion which is observed and (Do) is related 
to the amount of decoherence on the quantum system. The expecta- 
tion value of the back-reaction matrix (DP) quantifies the amount of 
back-reaction on the classical system. In the trivial case pe = 0, Eq. (20) 
places little restriction on the diffusion and Lindbladian rates appear- 
ing on the left-hand side. We already knew from refs. 28,29 that the pee 
must be a positive semi-definite matrix, and we also know that diffu- 
sion coefficients must be positive semi-definite. However, in the non- 
trivial case, the larger the back-reaction exerted by the quantum sys- 
tem, the stronger the trade-off between the diffusion coefficients and 
Lindbladian coupling. Equation (20) gives a general trade-off between 
observed diffusion and Lindbladian rates, but we can also find a trade- 
off in terms of a theory’s coupling coefficients alone. We show in the 
Appendix section “General trade-off between decoherence and diffu- 
sion coefficients” that the general matrix trade-off 


DTD" <2; (23) 


holds for the matrix whose elements are the couplings Dy ae 
for any CQ dynamics. Moreover, (I — DypDo pr = = 0, which tells us that 
Do cannot vanish if there is non-zero back-reaction. Equation (23) 
quantifies the required amount of decoherence and diffusion in order 
for the dynamics to be completely positive. In Eq. (23), and 
throughout, Do is the generalised inverse of pe, since pee is only 
required to be positive semi-definite. In the special case of a single 
Lindblad operator a =1 and classical degree of freedom, and when the 
only non-zero couplings are Dg := Do, D°,, := 2D, and D? , =1 this 
trade-off reduces to the condition D>Dp = 1 used in ref. 22. 

As a more general example, let us consider the class of theories 
that are continuous in phase space, and whose back- -reaction is gen- 


erated by a classical-quantum Hamiltonian A” which is only a 
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function of the canonical coordinates g;°. These are given by 


Bat Hes) He) 


bel pace! 1p on" Ol (24) 
2 gos} 0, if 0q; Q 0q; 


where H, is the purely classical Hamiltonian, p’ are the conjugate 
momenta, and Do and D; are q; dependent matrices with elements Do i 
and D>. Then the trade-off (23) implies that they must obey the matrix 
equation 2D,Dp <1. 

It is also useful to try to obtain an observational trade-off in terms 
of the total drift due to back-reaction as calculated in Eq. (11) 


> / dzTr/. Li Ly Q(Z)L a (25) 


pv#00 


It follows directly from Eq. (20) that when the back-reaction is 
sourced by either pe or ae we can arrive at the observational trade- 
off in terms of the total drift 


0 < 8(D,)(Dp) — (D})(DT)", Valz), (26) 
where the quantities appearing in Eq. (26) are now all observational 
quantities, related to drift, decoherence and diffusion as outlined in 
the previous subsection “Physical interpretation of the moments”. We 
believe that Eq. (26) should hold more generally, though we don’t have 
a general proof. 

In the case where the back-reaction is Hamiltonian at first order in 


the sense of Eq. (12), then Eq. (26) can be written as 


(0. SE) (aS) < 8(D;)(Do), Volz). (27) 


As a result, we can derive a trade-off between diffusion and 
decoherence for any theory that reproduces this classical limit and 
treats one of the systems classically. 

To summarise, whenever the back-reaction of the quantum sys- 
tem on the classical system induces a force on the phase space, then we 
have a trade-off between the amount of diffusion on the classical 
system and the strength of decoherence on the quantum system (or 
more precisely the strength of the Lindbladian couplings pi), This can 
be expressed both as a condition on the matrix of coupling co- 
efficients in the master equation, via Eq. (23) or in terms of observable 
quantities using Eqs. (20) and (26). In the case when the back-reaction 
is Hamiltonian, we further have Equation (27). We would like to apply 
this trade-off to the case of gravity in the non-relativistic, Newtonian 
limit. In order to do so, we will need to generalise the trade-off to the 
case of quantum fields interacting with classical ones, which we do in 
the subsection “Trade-off in the presence of fields”. The goal will be to 
understand the implications of treating the metric (or Newtonian 
potential) as being classical by using the trade-off when the quantum 
back-reaction induces a force on the gravitational field which, on 
expectation, is the same as the weak field limit of General Relativity. 


Trade-off in the presence of fields 

We would like to explore the trade-off in the gravitational setting and 
explore the consequences of treating the gravitational field as being 
classical and matter quantum. Since gravity is a field theory, we must 
first discuss classical-quantum master equations in the presence of 
fields. In the field-theoretic case, both the Lindblad operators and the 
phase space degrees of freedom can have spatial dependence, 
z(x),L,(x) and a general bounded CP map which preserves the 


classicality of the two systems can be written” 
p(z, th= : dz'dxdyA’" (z|z', t;x, VL OX, Z,2Z')oZ’, O)Li(y,z,2'), (28) 


where, as is usually the case with fields, in Eq. (28) it should be impli- 
citly understood that a smearing procedure has been implemented. 
We elaborate on some of the details when fields are introduced in the 
Appendix section “Classical-quantum dynamics with fields”. The con- 
dition for Eq. (28) to be completely positive on all CQ states is that for 
all vectors at x with components A,,(y, Z, 2’) 


/ dzdxdyA/,(x,z,2')M" (z|z';x, AZ, 2’) 20 (29) 


meaning that A“"(x, y) can be viewed as a positive matrix in pv and a 
positive kernel in x, y. In the field-theoretic case, one can still perform 
a Kramers-Moyal expansion and find a trade-off between the 
coefficients Do(x, y), D(x, y), D(x, y) appearing in the master equation. 
The coefficients now have an x,y dependence, due to the spatial 
dependence of the Lindblad operators. The coefficients 
D,(x, y), D2(x, y) still have a natural interpretation as measuring the 
amount of force (drift) and diffusion, whilst Do(x,y) describes the 
purely quantum evolution on the system and can be related to 
decoherence. 

Using the positivity condition in Eq. (29) we find the same trade of 
between coupling constants in Eq. (23) but where now D,(x, y) is the 
(p+ 1)n x (p + 1)n matrix-kernel with elements Do AX, y), D?"(x, y) is the 
(p + 1)n x p matrix-kernel with rows labelled by i, columns labelled by 
B, and elements bie ;(x,y), and Do(x, y) is the p x p decoherence matrix- 
kernel with ciements D(x, y). Here i€ {1,...,n}a€ {1,...,p} and 
ye {l,...,p+1}. In the field-theoretic trade off we are treating the 
objects in Eq. (23) as matrix-kernels, so that for any position- 
dependent vector b! (X), (Dob) (x) = f dy" y)b/(y), whilst for any 
position- dependent. vector (x), (Doa)*(x) = f dyp%8 (x, y)ag(). 
Explicitly, we find that positivity of the dynamics is equivalent to the 
matrix condition 


(30) 


br 
/ cay en ay | 2 D; a kA wT 


Dixy) Do(x,y) | a0) 


which should be positive for any position-dependent vectors bi,(X) and 
a,(x). This is equivalent to trade-off between coupling constants in 
Eq. (23) if we view (23) as a matrix-kernel equation. 

In order for the theory to be diffeomorphism invariant, we expect 
Do(x, y) and D2(x, y) to approach delta functions. We will not assume 
this, but we shall assume that the drift back-reaction is local, so that 
DP" (x, y) = 6(x, y)D°"(x). As we shall see in the next section, this is a 
natural assumption if we want to have back-reaction which is given by a 
local Hamiltonian. However, one might not want to assume that the 
form of the Hamiltonian remains unchanged to arbitrarily small dis- 
tances. With this locality assumption, Eq. (30) gives rise to the same 
trade-off of Eq. (23), where the trade-off is to be interpreted as a matrix 
kernel inequality. Writing this out explicitly we have 


/ dxedyal (x) DI4()(Do")aglX YD Y(X)ALX) S / dxdy2aj,(x)D5' (x, ya), 
(31) 


where asking that this inequality holds for all vectors a m(X) is equiva- 
lent to the matrix-kernel trade-off condition of Eq. (23). 

We give two examples of master equations satisfying the coupling 
constant trade-off in the Appendix section “Examples of Kernels 
saturating the decoherence diffusion coupling constants trade-off”. 
The decoherence-diffusion trade-off tells us how much diffusion and 
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stochasticity are required to maintain coherence when the quantum 
system back-reacts on the classical one. If the interaction between the 
classical and quantum degrees of freedom is dictated by unbounded 
operators, such as the mass density, then there can exist states for 
which the back-reaction can be made arbitrarily large. This is the case 
for a quantum particle interacting with its Newtonian potential 
through its mass density at arbitrarily short distances. Hence, if one 
considers a particle in a superposition of two peaked mass densities, 
then there can be an arbitrarily large response in the Newtonian 
potential around those points, and either there must be an arbitrary 
amount of diffusion, or the decoherence must occur arbitrarily fast. 
The former is unphysical, while the latter turns out to be the case in 
simple examples of theories such as those discussed in the Methods 
subsection “Decoherence rates”. 

Since our goal is to experimentally constrain classical-quantum 
theories of gravity, we shall hereby ask that the map (28) is CP when 
acting on all physical states p. If one allows for arbitrarily peaked mass 
distributions then the coupling constant trade-off of Eq. (31) should be 
satisfied. In the field-theoretic case, we can similarly find an observa- 
tional trade-off, relating the expected value of the diffusion matrix 
(D(x, y)) to the expected value of the drift in a physical state @ as we 
did in subsection “A trade-off between decoherence and diffusion’. 
This is done explicitly in the “Methods” subsection “Classical-quantum 
dynamics with fields’, using a field-theoretic version of the 
Cauchy-Schwartz inequality given by Eq.n (73), we find 


2(D,(x,x)) / dx’dy’ (Do(x',y)) = (D009) (De9)", G2) 


where Eq. (34) is to be understood as a matrix inequality with entries 


Do(x.y))= [ dar [DFPL—(%L}0)) (33a) 
(DP"(x,y)) = i dzTr [Dist yeool 4a], (33b) 
(D2Cy)\y= f d2Te [DY LaeL,O)]: (33c) 


Similarly, when the back-reaction is sourced by either po or pw it 
follows from Eq. (32) we can arrive at the observational trade-off in 
terms of the total drift due to back-reaction 


8(D9(x,X)) / dx'dy (Do (x,y) = (DEa))(DTe0)", G4) 


where 


(Oho), = | dzTe [DY paleo + DEL eC) + DM COL AL} 
(35) 


We shall now use the trade-off to study the consequences of 
treating the gravitational field classically. We will consider the back- 
reaction of the mass on the gravitational field to be governed by the 
Newtonian interaction (or more accurately, a weak field limit of Gen- 
eral Relativity). We shall then find that experimental bounds on 
coherence lifetimes for particles in superposition require large diffu- 
sion in the gravitational field in order to be maintained and this can be 
upper bounded by gravitational experiments. 

To summarise this section, we have derived the trade-off between 
decoherence and diffusion for classical-quantum field theories, both 
in terms of coupling constants of the theory and in terms of 


observational quantities. This trade-off puts tight observational con- 
straints on classical theories of gravity which we now discuss. 


Physical constraints on the classicality of gravity 

In this section, we apply the trade-off of Eq. (30) to the case of gravity. 
A number of classical-quantum models of Newtonian gravity have 
been proposed”**’, but since the trade-offs derived in the previous 
section depend only on the back-reaction, or drift term, they are 
insensitive to the particulars of the theory. We shall consider the 
Newtonian, non-relativistic limit of a classical gravitational field which 
we reproduce in the “Methods” subsection “Newtonian limit of CQ 
theory”. A fuller discussion, including a derivation of the Newtonian 
limit starting from the covariant theories of refs. 11,39 can be found in 
ref, 81. It is in taking this limit where some care should be taken, since 
one is gauge fixing the full general relativistic theory. We denote ® to 
be the Newtonian potential and in the weak field limit of General 
Relativity, it has a conjugate momenta we denote by 77m. We assume: 
(i) The theory satisfies the assumptions used to derive the master 
equation as in subsections “Trade-off in the presence of fields”; in 
particular that the theory be a completely positive norm- 
preserving Markovian map, and that we can perform a short- 
time Kramers-Moyal expansion as in “Methods” subsection 
“Classical-quantum dynamics with fields”. 

We apply the theory to the weak field limit of General Relativity, 
whereas recalled in “Methods” subsection “Newtonian limit of CQ 
theory” the Newtonian potential interacts with matter through its 
mass density m(x), 


(ii) 


H,(®)= / PxOx)m(x). (36) 
and the conjugate momentum to © satisfies 
Tg = ve — m(x), (37) 


where in the c > ~ limit the momentum constraint 77m ~ 0 is imposed 
and we recover Poisson’s equation for the Newtonian potential. We 
assume this limit of General Relativity is satisfied on expectation, at 
least to leading order. This may be an overly strong assumption, since 
the weak field limit may cease to be valid at short distances when the 
diffusion becomes large. A relativistic treatment is initiated in (J. 
Oppenheim and A. Russo, manuscript in preparation). It is also worth 
noting that General Relativity has not been tested at distances shorter 
than the millimeter scale, and here we assume it holds to arbitrarily 
short distances. 
(iii) In relating Dg to the decoherence rate of a particle in super- 
position, we shall assume that the state of interest is well 
approximated by a state living in a Hilbert space of fixed particle 
number. We believe this is a mild assumption: ordinary non- 
relativistic quantum mechanics is described via a single particle 
Hilbert space, and we frequently place composite massive 
particles in superposition and they do not typically decay into 
multiple particles. 
We will assume that the diffusion kernel D,(®,x,x’) does not 
depend on 79 i.e. it is minimally coupled. This is reasonable, since 
in the purely classical case matter couples to the Newtonian 
potential and not its conjugate momenta. 


(iv) 


With these assumptions, and treating the matter density as a 
quantum operator m(x), this tells us that in order for the back-reaction 
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term to reproduce the Newtonian interaction on average 


Tr[{AH,,0}] =Tr I/ Bximoo ee | 


atin be; (38) 
== Trl [ xDE" (®, Teo, X)Ly(X)<—L}(x) |, 
2 LP, Ho, Xk yl li, (x) 
then we must pick 

(Dy, 1,(,Mo,X)) = — (h(X)), (39) 


meaning that the back-reaction matrix Di is nonvanishing. In the 
“Methods” subsection “Newtonian limit of CQ theory” we give 
examples of master equations for which Eq. (39) is satisfied, but their 
details are irrelevant since we only require the expectation of the back- 
reaction force to be the expectation value of the mass—a necessary 
condition for the theory to reproduce Newtonian gravity. 

As a consequence of the coupling constant and observational 
trade-offs derived in Eqs. (31) and (32), a non-zero D,, its implies that 
there must be diffusion in the momenta conjugate to 7. This diffusion 
is equivalent to adding a stochastic random process /(x, ¢) (the Lan- 
gevin picture), to the equation of motion (37) to give 


2 


bs ose (40 
"0 4G : 


= m(x) + u®, mit, x), 


where we allow some colouring to the noise via a function u(®, m) 
which can depend on 9, and the matter distribution m (assumption 
(iv)). The noise process satisfies 


mol, 0)=0, En olWe, Ow(y,t)|=2(D,,y,9)) 66,0), (Al) 
where we have defined (D,(x, y, ®)) = Tr[D5"(x,y, O)L,(X)PLi(y)], and p 
is the quantum state for the decohered mass density. Here the m, 
subscripts of E,,@ allow for the possibility that the statistics of the 
noise process can be dependent on the Newtonian potential and mass 
distribution of the particle. The restriction on E,, o[w/(x, t)] follows 
from assumption (ii). If u/(x, 0) is Gaussian, Eq. (41) completely deter- 
mines the noise process, but in general, higher-order correlations are 
possible, although they need not concern us here, since we are only 
interested in bounding the effects due to D2(x, y, ®). 

In the non-relativistic limit, where c >, we impose the momen- 
tum constraint 70 and we recover Poisson’s equation for gravity, 
but with a stochastic contribution to the mass. This is precisely as 
expected on purely physical grounds: in order to maintain coherence 
of any mass in superposition, there must be noise in the Newtonian 
potential and this must be such that we cannot tell which element of 
the superposition the particle will be in, meaning the Newtonian 
potential should look like it is being sourced in part by a random mass 
distribution. In other words, the trade-off requires that the stochastic 
component of the coupling obscures the amount of mass m at the 
different points in space where the mass may be found. 

In the case where u is independent of ®, it is simple to solve 
Eq. (40) in terms of Green’s function for Poisson’s equation as in the 
“Methods” subsection “Detecting gravitational diffusion’. A formal 
treatment of solutions to non-linear stochastic integrals of the more 
general form of Eq. (40) can be found in ref. 82. A higher precision 
calculation would involve a full simulation of CQ dynamics, for 
example using unravelling methods”*”’ or the path integral as in ref. 81. 
However, care should be taken, as we have found that relativistic 
corrections put constraints on the degree of diffusion even at low 
energy (J. Oppenheim and A. Russo, manuscript in preparation), and 
one should bear this in mind when drawing conclusions on the models 
presented here. 


In*° it was shown that there are two classes of CQ dynamics, at 
least in the sense that there are those with continuous trajectories in 
phase space and those which contain discrete jumps. For the class of 
continuous CQ models (see ref. 24 and Appendix section “Continuous 
master equation”), we know that J(x, t) should be described by a white 
noise process in time, and its statistics should be independent of the 
mass density of the particle. 

For the discrete class (see ref. 11 and J. Oppenheim, “The con- 
straints of a continuous realisation of post-quantum-classical gravity", 
manuscript in preparation) and “Methods” subsection “Discrete mas- 
ter equation’), /(x, t) can involve higher order moments, and will gen- 
erally be described by a jump process”®”’. Its statistics can also depend 
on the mass density, since in general the diffusion matrix Dey couples 
to Lindblad operators. It is worth noting that the discrete CQ theories 
considered in"”°*’ generically suppress higher order moments, and 
often we expect that we can approximate the dynamics by a Gaussian 
process, but this need not be the case in general. 

The stochastic contribution to the Newtonian potential leads to 
observational consequences which can be used to experimentally test 
and constrain CQ theories of gravity for various choices of kernels 
appearing in the CQ master equation. One immediate consequence is 
that the variation in Newtonian potential leads to a variation of force 
experienced by a_ particle or composite mass _ via 
Foz —f d?xm(x)V@(x). We can also estimate the time-averaged 
force via Ot a F tot Where AT is the time resolution over which the 
force is measured and is the useful quantity when compared with 
experiments. In the “Methods” subsection “Table-top experiments” we 
impose the constraint 7 = O in Eq. (40) and find that the variance of the 
magnitude of the time-averaged force experienced by a particle in a 
Newtonian potential is given by Eq. (146), 


2G 8 AS a8 48a 

of= Fp | Pxdyd x’ dP’ y'm(x)m(y) 

ts ts 42 
(¥-*X)-G-7V) 
Ix-x'Byy-yP 


(Dx(x',Y,®)), 


where the variation is averaged over the time resolution AT. We will use 
this to estimate the variation in precision measurements of mass, such 
as modern versions of the Cavendish experiment for various choices 
of (D,(x’, y’, ®)). 

On the other hand, experimentally measured decoherence rates 
can be related to Do. The important point is that the decoherence rate 
is dominated by the background Newtonian potential ®, due to the 
Earth. In the “Methods” subsection “Decoherence rates”, we show that 
for a mass whose quantum state is a superposition of two states |L) and 
|R) of approximately orthogonal mass densities m (x), mp(x), and 
whose separation we take to be larger than the correlation range of 
Do(x, y), the decoherence rate is given by 


1 + F 
= 5 | dedyo efx, WALILjOMLACOIL) + (RLZONLAOIR)). (43) 


Via the coupling constant trade-off, Eqs. (42) and (43) then give 
rise to a double-sided squeeze on the coupling D>. Equation (42) upper 
bounds D, in terms of the uncertainty of acceleration measurements 
seen in gravitational torsion experiments, whilst the coupling constant 
trade-off Eq. (43) lower bounds D, in terms of experimentally mea- 
sured decoherence rates arising from interferometry experiments. 

We now show this for various choices of diffusion kernel, with the 
details given in the “Methods” subsection “Table-top experiments”. 
The bounds are summarised in Table 1. The diffusion coupling strength 
will be characterised by the coupling constant D2, which we take to bea 
dimension-full quantity with units kg” sm°, and is related to the rate of 
diffusion for the conjugate momenta of the Newtonian potential. We 
upper bound D, by considering the variation of the time averaged 
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Table 1 | Current experimental bounds on classical-quantum theories for different master equations and functional depen- 


dence on the diffusion coefficient 


Master equation Diffusion kernel 


Experimental squeeze 


Continuous (ultra-local, non-rel.) 
Continuous (Eq. (116) or (118)) 


DO; x, y) = Da()6(x, y)D2(®) = d,c°" 
Da(®;x,y)= — D(@)V76(x,y)Do() = 7,07" 


10“ > D, > 10° kg? sm® (Eq. (44)) 
10° > 3D, = 10-* kg’?sm- (Eq. (47) 


Discrete (ultra-local) 


Do(@; x,y) = Fe Do()4(X,y)D2(@) =>>,07 0” 


10-1 > #2 > 10° kg (Eq. (46)) 


The diffusion coefficient is bounded from above by observed acceleration variations 02 seen in precision mass experiments via Eq. (42). In all cases the master equation is assumed to saturate the 
bound which is used to find the lower bound the amount of diffusion on the quantum system by bounding Do from coherence rates via Eq. (43). It is seen that minimally coupled continuous models 
which are non-relativistic and do not create spatial correlations (we call these ultra-local) and have polynomial dependence on the Newtonian potential are ruled out, while continuous models with 
non-local correlations, such as the Diosi-Penrose (DP) kernel of Eq. (116)or the relativity inspired kernel of Eq. (118), and ultra-local discrete models are less constrained. Here lp, mp are the Planck 
length and Planck mass respectively, which are required in order for the dimensions of D2(®) to be the same in all cases. 


acceleration o,, = $5 for a composite mass M which contains N atoms 
which we treat as spheres of constant density p with radius ry and mass 
my. We lower bound D) via the coupling constant trade-off of Eq. (30) 
and then by considering bounds on the coherence time for composite 
particles with total mass M, and which are made up of N, constituents, 
each with typical length scale when in superposition R, and volume Vy. 

For continuous dynamics (D3(x, y, ®)) = D2(x, y, D) since the dif- 
fusion is not associated with any Lindblad operators. Let us now con- 
sider a very natural kernel, namely D(x, y; ®) = D2(®)6(x, y) which is 
both translation invariant, and does not create any correlations over 
space-like separated regions. We call dynamics which does not create 
correlations over space-like separated regions ultra-local since the- 
ories that are not of this form can still be non-signalling. This is a 
natural kernel from the point of view of constructing theories which 
are diffeomorphism invariant. We also label this model as being non- 
relativistic since it does not include various relativistic corrections to 
the diffusion. 

The decoherence rate for this kernel is found in the “Method” 
subsection “Decoherence rates” and follows immediately from Eq. 
(137). For a nucleon of mass M, and wavepacket volume V4, it is 
A=2D)M3/V,. In general, the squeeze will depend on the functional 
choice of D2(®) on the Newtonian potential. However, in the presence 
of a large background potential ®,, such as that of the Earth’s, we will 
often be able to approximate D2(®) = D2(®,). This is true for kernels 
that depend on ® and V9, though the approximation does not hold 
for all kernels, for example D2 ~—-V?® of Eq. (118) which creates diffu- 
sion only where there is mass density. For diffusion kernels D(®,) 
where the background potential is dominant, we find the promised 
squeeze on D2(®p) 
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where V, is the volume of space over which the background Newtonian 
potential is significant. V, enters since the variation in acceleration is 
found to be 


2 D.C 


. d?x’D,(®,), 
°a~ nat | 9% 22o) 


(45) 


where the d°x’ integral is over all space. 

This immediately rules out continuous theories of Newtonian 
gravity with noise everywhere, i.e., with a diffusion coefficient inde- 
pendent of the Newtonian potential, since the integral will diverge. We 
consider the relativistic case elsewhere. 

Standard Cavendish-type classical torsion balance experiments‘? 
measure accelerations of the order 10’ ms” over minutes AT - 107, so 
a very conservative bound is 0, ~ 107’ ms~, whilst for a kg mass N ~ 107° 
and ry~10"'m. Conservatively taking V,~1r2hm> where rg is the 
radius of the Earth and hf is the atmospheric height gives 
D,<10“ kg’? sm™*. The decoherence rate A is bounded by various 
experiments®. Typically, the goal of such experiments is to witness 
interference patterns of molecules that are as massive as possible. 


Taking a conservative bound on A, for example, that arising from the 
interferometry experiment of which saw coherence in large organic 
fullerene molecules with total mass M,=10™ kg over a timescale of 
0.1s, gives an upper bound on the decoherence rate A<10's™. Full- 
erene molecules are made up of Nj ~ 10° particles with a typical atomic 
size 10° m. After passing through the slits the molecule becomes 
delocalised in the transverse direction on the order of 10 m before 
being detected. Since the interference effects are due to the super- 
position in the transverse x direction, which is the direction of align- 
ment of the gratings, it seems like a reasonable assumption to take the 
size of the wavepacket in the remaining y, z direction to be the size of 
the fullerene, since we could imagine measuring the y,z directions 
without effecting the coherence. We, therefore, take the volume 
V,~1010%107 m? =10’ m?, which gives D> >10~° kg? sm°, and sug- 
gests that classical-quantum theories of Newtonian gravity with ultra- 
local continuous noise are ruled out by experiment. 

On the other hand, the discrete models appear less constrained 
due to the suppression of the noise away from the mass density. For 
example consider the ultra-local discrete jumping models, such as the 
one given in the section “Discrete master equation” which have 


(Do(x,y, Dy)) = HPs) mx), where mp = YE: is the Planck mass and 


[p= ys is the Planck length, required to ensure D, has the units of 
kg’ sm. We find the squeeze on D5 


OANrVAT . Ip p,> Ma 


> . (46 
MyG? Mp A 


and plugging in the numbers tells us that discrete theories of classical 
gravity are .not ruled out by experiments and we 
find 10" kg> ;£-D, >10-* kg. 

We can also consider other noise kernels, with examples and some 
discussion is given in the section “Examples of Kernels saturating the 
decoherence diffusion coupling constants trade-off’. A natural kernel 
is D,(x,y, D,)= — [2D(®,)V76(x, y). The inverse Lindbladian kernel 
satisfying the coupling constants trade-off is to zeroeth order in M(x), 
the Diosi-Penrose kernel D(x, y, Oy) = Poe We here consider higher- 
order terms such as those coming from the relativistic theory and in 
particular the diffusion kernel of Eq. (118). For this choice of dynamics, 
we find the squeeze for D2 in terms of the variation in acceleration 


o2Nr3, AT 


oh 2 GD32 (47) 


Using the same numbers as for the ultra-local continuous model, 
with R,~V}/?~10-" m we find that classical torsion experiments 
upper bound D, by 10~° kg’sm-! > [2D>, whilst interferometry 
experiments bound D, from below via iD, >10-* kg” sm-}. 

Equations (44), (46) and (47) show that classical theories of gravity 
are squeezed by experiments from both ways. We have here been 
extremely conservative, and we anticipate that further analysis, as well 
as near-term experiments, can tighten the bounds by orders of 
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magnitude. There are several proposals for table-top experiments to 
precisely measure gravity, some of which have recently been per- 
formed, and which could give rise to tighter upper bounds on D>. Some 
of these experiments involve millimeter-sized masses whose gravita- 
tional coupling is measured via torsional pendula®’’, or rotating 
attractors’. With such devices, the gravitational coupling between 
small masses can be measured while limiting the amount of other 
sources of noise. There are proposals for further mitigating the noise 
due to the environment, including inertial noise, gas particle collisions, 
photon scattering on the masses, and curvature fluctuations due to 
other sources***°, Other experiments are based on interference 
between masses; for example, atomic interferometers allow for the 
measurement of the curvature of space-time over a macroscopic 
superposition®’ *°. 

We can get stronger lower bounds via improved coherence 
experiments. Typically, the goal of such experiments is to witness 
interference patterns of molecules that are as massive as possible, 
while here, we see that the experimental bound on CQ theories is 
generically obtained by maximising the coherence time for massive 
particles with as small wave-packet size V,. 

Thus far we have considered local effects on particles due to 
diffusion. While this enables us to rule out some types of theories, the 
bounds are generally weak if one wants to rule out all of them. How- 
ever, it may be possible to do so via cosmological considerations. In 
attempting to place experimental constraints on this diffusion, it is 
also worth considering other regimes, such as longer range effects 
which might be detected by gravitational wave detectors such as LIGO, 
or table-top interferometers’””". We leave a detailed study of the effect 
of gravitational diffusion on cosmological scales and LIGO to future 
work. It suffices to mention that the effect will again depend on the 
form of the kernel D,(x, x’). Our estimates (J. Oppenheim and Z. Weller- 
Davies, “Estimating space-time diffusion in interferometers", unpub- 
lished note) suggest that local effects from table-top experiments 
currently place a stronger bound on gravitational theories than LIGO 
currently does. In particular, unlike gravitational wave measurements, 
which are reasonably high-frequency events requiring extraordinarily 
high precision in the relative displacement of the arm length from its 
average, it is preferential to have a lower precision measurement, 
which occurs over a longer time period to allows for the diffusion in 
path length to build up, and with a smaller uncertainty in the average 
length of the arm itself. Furthermore, since the LIGO arm is kept in a 
vacuum, we do not expect strong bounds on discrete models where 
the diffusion is associated with an energy density. 


Discussion 
Anumber of direct proposals to test the quantum nature of gravity are 
expected to come online in the next decade or two. These are based on 
the detection of entanglement between mesoscopic masses inside 
matter-wave interferometers” ”°’””. For these experiments, some 
theoretical assumptions are needed: one requires that it is only grav- 
itons that travel between the two masses and mediate the creation of 
entanglement. If this is the case, then the onset of entanglement 
implies that gravity is not a classical field. These can be thought of as 
experiments that if successful, would confirm the quantum nature of 
gravity (although other alternatives to quantum theory are possible”’). 
Here, we come from the other direction, by supposing that gravity 
is instead classical, and then exploring the consequences. Theories in 
which gravity is fundamentally classical were thought to have been 
ruled out by various no-go theorems and conceptual difficulties. 
However, these no-go theorems are avoided if one allows for non- 
deterministic coupling as in"?"?°°°?”, We have here proven that this 
feature is indeed necessary and made it quantitative by exploring the 
consequences of complete positivity on any dynamics that couples 
quantum and classical degrees of freedom. Complete positivity is 
required to ensure the probabilities of measurement outcomes remain 


positive throughout the dynamics. We have shown that any theory 
which preserves probabilities and treats one system classically is 
required to have fundamental decoherence of the quantum system, 
and diffusion in phase space, both of which are signatures of infor- 
mation loss. Using a CQ version of the Kramers—Moyal expansion, we 
have derived a trade-off between decoherence on the quantum sys- 
tem, and the system’s diffusion in phase space. The trade-off is 
expressed in terms of the strength of the back-reaction of the quantum 
system on the classical one. We have derived the trade-off both in 
terms of coupling constants of the theory and in terms of observa- 
tional quantities that can be measured experimentally. 

In the case of gravity, the trade-off places a lower bound on the 
rate of diffusion of the gravitational degrees of freedom in terms of the 
decoherence rate of particles in superposition. We find that theories 
that treat gravity as fundamentally classical, are not ruled out by cur- 
rent experiments, however, we have been able to rule out a broad 
parameter space of Newtonian theories. This is done partly through 
table-top observations via Eqs. (44), (46) and (47). Given any diffusion 
kernel, we can compute the inaccuracy of mass measurements due to 
fluctuations in the gravitational field, and using the trade-off, we can 
derive a bound on the associated decoherence rate. This allows us to 
rule out broad classes of theories in terms of their diffusion kernel. For 
example, we are able to rule out a number of non-relativistic theories 
which back-react continuously in phase space. 

Any theory that treats gravity classically has fairly limited freedom 
to evade the effects of the trade-off. There is the freedom to choose the 
diffusion or decoherence kernels D,(x, x’) and Do(x,x’), but the trade- 
off restricts one in terms of the other. Then, because of the results 
proven in®°, one can consider two classes of theory, those which are 
continuous realisations and whose diffusion can only depend on the 
gravitational degrees of freedom, and discrete theories whose diffu- 
sion can also depend directly on the matter fields. Examples of both 
classes of the theory are given in the “Methods” subsection “Newtonian 
limit of CQ theory”. Finally, one could consider theories that do not 
reproduce the weak field limit of General Relativity to all distances, 
namely we could imagine that the interaction Hamiltonian of Eq. (36) 
does not hold to arbitrarily short distances, or arbitrarily high mass 
densities. This is reasonable since we do expect the Newtonian theory 
to break down at short distances where relativistic corrections at high 
energy affect the low-energy behaviour of the theory. One could also 
consider modifying D,(x,x’) in some other way, for example, by 
making it slightly non-local, or by disallowing arbitrarily high mass 
densities, or by including an additional contribution such as the fric- 
tion term discussed in the continuous master equation. All of these 
modifications would seem to violate Lorentz invariance in some way, 
and likely lead to observational consequences”’. 

Here, we have only given an order of magnitude estimate of when 
gravitational diffusion will lead to appreciable deviations from New- 
tonian gravity. The most promising experiments bounding the diffu- 
sion appear to be table-top experiments which precisely measure the 
mass of an object. This is an area that is important from the perspective 
of weight standards, for example, those undertaken by NIST on the 1 kg 
mass standard K20 and K4”*. The increased precision and measuring 
time of Kibble Balances® and atomic interferometers*”**°°”” would 
make such measurements an ideal testing ground, both to further 
constrain the diffusion kernel and to look for diffusion effects, whose 
dependence on the test mass is outlined in the “Methods” subsection 
“Detecting gravitational diffusion’. Here, we have found that the 
resolution time AT over which variations of acceleration are estimated 
affects the strength of the bound, and it would be helpful if future 
experiments reported this value. Since we have found that CQ theories 
predict an uncertainty in mass measurements it is perhaps intriguing 
that different experiments to measure Newton’s constant G yield 
results whose relative uncertainty differs by as much as 
5x10* m’?kg‘s”, which is more than an order of magnitude larger 
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than the average reported uncertainty” **. If one were to try and 
explain the discrepancy in G measurements via gravitational diffusion, 
then for all the kernels we studied in Section ’Physical constraints on 
the classicality of gravity’ we find that the variation in acceleration 
depends on — the number of nucleons in the test mass, so that masses 
with smaller volume should yield larger uncertainty and this would be 
the effect to look for in measurement discrepancies. The relatively 
large uncertainty in such measurements, also makes it challenging for 
table-top experiments to place strong upper bounds on gravitational 
diffusion. 

Turning to the other side of the trade-off, improved decoherence 
times would further squeeze theories in which gravity remains classi- 
cal. While a current experimental challenge is to demonstrate inter- 
ference patterns using larger and larger mass particles, we here find 
the bounds in Eqs. (44) and (46) depend on the expectation of the 
particle’s mass density in ways that depend on the particular kernel. 
Thus interference experiments with particles of high mass density 
rather than mass can be preferable. There are also kernels, for which 
the relevant quantity is the expectation of the mass density or M7/V,, 
which will depend on both the particle’s mass M, and volume V, of the 
wave-packet used in the interference experiment, a quantity which is 
not always obtainable from many reports on such experiments. While 
this dependence might initially appear counter-intuitive, it follows 
from the fact that in order to relate the trade-off in terms of coupling 
constants to observational quantities, and in particular, the deco- 
herence rate, we took expectation values of the relevant quantities to 
get a trade-off in terms of only averages. And indeed the decoherence 
rate, which is an expectation value, can easily depend on the wave- 
packet density, as we see from examples in the section “Deco- 
herence rates”. 

Since we here show that all theories that treat gravity classically 
necessarily decohere the quantum system, another constraint on 
theories that treat gravity classically is given by constraints on funda- 
mental decoherence. These are usually constrained by bounds on 
anomalous heating of the quantum system’®. However, these con- 
straints are not in themselves very strong, since fundamental deco- 
herence effects can be made arbitrarily weak. In the simplified model in 
the “Methods” subsection “Newtonian limit of CQ theory”, the strength 
of the decoherence depends on the strength of the gravitational field, 
thus, constraints due to heating’*™ can be suppressed, either by 
scaling the Lindbladian coupling constants or by having strong deco- 
herence effects more pronounced near stronger gravitational fields 
such as near black holes where one expects information loss to occur. 
The necessity for decoherence to heat the quantum system is further 
weakened by the fact that the dynamics are not Markovian on the 
quantum fields, if one integrates out the classical degrees of freedom, 
space-time acts as a memory. This potentially captures some of the 
non-Markovian features advocated in ref. 112, who recognised that 
Markovianity is a key assumption in attempts to rule out fundamental 
decoherence or information loss. Here, however, we see that there is 
less freedom than one might imagine. If the Lindbladian coupling 
constants are made small to reduce direct heating, the gravitational 
diffusion must be large. Thus, heating constraints which place bounds 
on Do(x,x’) place additional constraints on D,(x,x’). In®, it was found 
that for the Newtonian models of ref. 33, large D,(x,x’) creates sec- 
ondary heating which further constrain the theory experimentally. The 
decoherence-diffusion trade-off implies that this is a general feature of 
all theories which treat gravity classically. 

While the absence of diffusion could rule out theories where 
gravity is fundamentally classical, the presence of such deviations, at 
least on short time scales, might not by itself be a confirmation of the 
classical nature of gravity. Such effects could instead be caused by 
quantum theories of gravity whose classical limit is effectively described 


by Oppenheim”. In other words, one might expect some gravitational 
diffusion, because, from an effective theory point of view, one is in a 
regime where space-time is behaving classically. There are even claims 
that holographic effects could cause stochasticity’* "> in the gravita- 
tional field. However, the trade-off we have derived is a direct con- 
sequence of treating the background space-time as fundamentally 
classical. In a fully quantum theory of gravity, the interaction of the 
gravitational field with particles in a superposition of two trajectories 
will cause decoherence, but coherence can then be restored when the 
two trajectories converge. This is because the particle’s position is 
entangled with the gravitational field (or dressed by it), and this 
entanglement is erased when the different paths of the superposition 
converge. This is what happens when electrons interact with the elec- 
tromagnetic field while passing through a diffraction grating, yet still 
form an interference pattern at the screen. This is a non-Markovian 
effect—the which-path superposition decoheres almost immediately, 
but this is false-decoherence'’ so the amount of diffusion can be arbi- 
trarily small and is unrelated to the coherence time of the superposition. 

On the other hand, the trade-off we derived is a direct con- 
sequence of the positivity condition, which is a direct consequence of 
the Markovian assumption. In the non-Markovian theory where Gen- 
eral Relativity is treated classically, one still expects the master equa- 
tion to take the form found in”, but without the matrix whose elements 
are D*” needing to be positive semi-definite at all times'’”""®. It would 
therefore be surprising, if a quantum theory of gravity predicted 
anything close to the level of diffusion predicted by the decoherence- 
vs-diffusion trade-off, as there would be no need for diffusion to 
explain the coherence of superpositions. The regime in which the 
classical-quantum theory can be regarded as an effective one is taken 
up in ref. 119, both to address the issue of false decoherence, and also 
to explore the regime in which the classical-quantum theory may be a 
useful tool to understand the back-reaction of quantum matter in 
space-time, such as during black-hole evaporation, and during infla- 
tion. If we instead regard the theory as describing a fundamentally 
classical space-time, then it follows from the decoherence-diffusion 
trade-off, that the diffusion is either fundamental or its source is not 
describable within quantum or classical mechanics J. Oppenheim, 
“Post-quantum soup", unpublished note). 


Methods 
Positivity conditions and the trade-off between decoherence 
and diffusion 
In this section, we will introduce two forms of positivity conditions 
used to prove the decoherence diffusion trade-off. 

The first inequality we would like to introduce is 


i dzA),(z,z’)A"" (z|z', 6t)A,(z,z')20, (48) 


which holds for any A,,z, z’) for which Eq. (48) is well defined: i.e., so 
that the distributional derivatives in Eq. (48) are well defined. 

We can derive the positivity condition (48) from the positivity of 
A“’(z|z'), which must be a positive semi-definite matrix in zv. More 
precisely, the eigenvalues of A“’(z|z’), which we denote by A“(z|z’) 
must be positive. They must be positive in the distributional sense, 
since we allow for the case that A“(z|z’) is a positive distribution, for 
example M(z\z’) ~ 6(z — z’). Hence we require 

/ dzdz‘A"(z\z')P(z,z’) (49) 


is positive for any positive smearing function P(z,z’). Since each A” 
must be positive, we can also pick a different smearing function for 
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each py, so that 


/ dzdz'M*(z\z')P,(z,2’) (50) 


should be positive for any vector P,,(z, 2’) with all positive entries. We 
can then write the matrix A’”(z|z’) in terms of its eigenvalues 


AY (z|2') = US" (z\z' AP (zz U5 (z12’). (51) 


We can then see the positivity of Eq. (48) directly since 


/ dzA),(z,z')AY(z|z', 6t)A,(z,2')= fi dz(UA)),(z|z')AK(z|z' (UA) (2, 2’) 


= / dz\(UA),|7(z, 2’ AM(z|z’) 
(52) 


which is positive as a consequence of Eq. (50). 
As a consequence of Eq. (48) being positive, we also know that 


Tr / dzM""(z|z')0,(z, z\pt2')0\(2,2') >0 (53) 


will be positive for any vector of operators (potentially phase space 
dependent) O,(z,z’). This follows from the cyclicity of the trace and 
the fact that NY (z|z')ON(z, 2')0,(Z,z') will be a positive operator so 
long as Eq. (48) holds. A common choice of O, would be the Lindblad 
operator L,, appearing in the master equation. 

The inequality in Eq. (48) proves useful to derive positivity 
conditions on the coupling constants appearing in the master 
equation, whilst Eq. (53) is useful in deriving the observational 
trade-off for the continuous master equation as we shall now 
discuss. 


General trade-off between decoherence and diffusion 
coefficients 
We can get a general trade-off between the decoherence and diffusion 
coefficients which appear in the master equation, arriving at a trade-off 
between the decoherence and diffusion coefficients in terms of the 
back-reaction drift coefficient Di“. 

Consider Eq. (48), and choose A, =6),dq + b,(Z — 2’); By inte- 
grating parts over the phase space degrees of freedom, we find 


201, DS” bi, + bi Diag + a, DiXbi, + a, DE" ag > 0. 


iv Op a 


Taking i€ {1,...,n}a € {1,..., p} and we {1,...,p +1}, we can write 
this as a matrix positivity condition 


» .[2D, D|{b 
(eal 0 (55) 
pp, |La 
where Dj is the (p+1)nx (p+ 1)n matrix with elements DY, DY" is 


the (p+1)nxp matrix with rows labelled by yi and columns 
labelled by 6 with elements pie and Do is the px p decoherence 
matrix with elements D??. fa describes the quantum back- 
reacting components of the drift. Equation (55) is equivalent to 
the condition that the ((p + 1)n+ p) x((p + 1)n+ p) matrix 


2 9 


(56) 
DE Do 


Since we know D3 and Do must be positive semi-definite, we know from 
Schur decomposition that 

2D, >D*"D, De", (57) 
and (I — DoDo — 1p" =0, where Do is the generalised inverse of Do. 
Furthermore, if Do vanishes, then clearly D°" must also vanish in order 
for (56) to be positive semi-definite. 


Classical-quantum dynamics with fields 

In this section, we describe CQ dynamics in the case where the Lind- 
blad operators and the phase-space degrees of freedom can have 
spatial dependence z(x), L,,(x). 

For the case of fields, operators O(x) constructed out of local 
fields @(x) will in general be unbounded and hence the Stinespring 
dilation theorem does not hold. This problem is a common one in the 
study of algebraic quantum field theory and we can get around it by 
considering the case in which operators of interest are obtained by 
smearing the local fields over bounded functionals F. For example, we 
can first smear the local field fields over a smearing function 
f, @¢=Jdx@(x)fo) and then consider bounded functions of @; such as 
F(@,) =e. In doing this we can write a CQ version of the Stinespring 
dilation theorem exactly and proceed along the lines of Oppenheim" 
to show that any completely positive CQ map can be written in the 
form 


p (z= i dzdxdyA"’(z|z':x, y)L 0% Z, Zaz (y, 2,2"), (58) 
where the positivity condition states 
/ dzdxdyA/,(x,Z,2')M"(z|z';x, Y)A,(,Z,2') 20. (59) 


We shall assume that we deal with dynamics which can be written 
in Lindblad form, as is usually assumed in the unbounded case”. 


CQ Kramers-Moyal expansion for fields 
Just as in the section “The CQ Kramers-Moyal expansion”, we can 
formally introduce the moments of the transition amplitude 


MY 


1, iy...5in 


(Wy, ... Wy: X,Y, 50) = / DzN" (z\z';x,y, 6t) (2 — 2"), (W,)...(Z —2');, (Wp) 
(60) 


which we assume to exist; which might involve a smearing of the 
operators z(x). Defining Lo(x)=6(x)I, we can define the coefficients 


Dr',,.i, implicitly via 
Mnni,i, (2 1 Wy ++ WysX,Y, Ot) = 655 + OtN'D yp (Wy, --- WyiX,Y, OC). 


(61) 


The characteristic function then takes the form 
cm uzix.y= [ Dre lammemsnnmea'sxy) (62) 


and expanding out the exponential this takes the form 


Vv /, _ a ¢ yw 
cH wzixn= > | dw,...dw, a M 


Nn, iy... 


i,(2Z Wy, Wy; X,y, Ot) 
(63) 


performing the inverse Fourier transform, allows us to write the 
transition amplitude in terms of functional derivatives of the delta 
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function 


ME cig(Z Wir «Wy X,Y, 0) o 


NM” (z|z';x,y, Ot) = [aw ... dw, ri 
GIz;x,y, 06) > : n n! 6Z;, (Wy) ...2), (Wp) 
(64) 
and we can use this to write a CQ master equation in the form 
doz, 6t)_ = Z 6" E 
a yy) dw, ...dw,(—1) Bea zt nt ‘aig ZW - -Wy)Q(2)) — i{H, a2] 


1 : 
+ / dxdyD™ (z; x, y)La(X)ol ZL p(y) — 5 DO (zx (Lge), @) 


+ yf dxayauy.. -dw,(—1)" 


n=0 p00 


Diy ig(Z Wr -- 


6" ; 
ACHR WpiX WlyOOLLO)). 
i, aad Fs 


(65) 


Since we are interested in studying dynamics with local back- 
reaction, we shall hereby take D{"(z,w;x,y)=Di" (~)6(x, G(x, w). 
By the decoherence diffusion trade-off, which we derive in the 
next subsection, this also means that the diffusion matrix 
Dy’ (Z,W ,W2,X,y) is lower bounded by the matrix 

Di (oo(D5! Jap, y)Dhv" (y)6(w,,x)6(W>,y). This can be seen more pre- 
geal: by taking Eq. (59) with A,(x)= On Ag(X) + f dub, (x, W)\(Z — 
z'\(x, W) and applying the same methods as ‘in the subsection “Trade- 
off between diffusion and decoherence couplings in the presence of 
fields’. Without loss of generality we thus _ take 
D2(Z, Wy, W2, X, Y) = Do(Z, xX, Y)6(x, Wy)6(y, W2). 


Trade-off between diffusion and decoherence couplings in the 
presence of fields 

In the field-theoretic case, the positivity condition is given by Eq. (59) 
and we can find a trade-off between decoherence and diffusion by 
considering A,,(x) = 571 Aq(X) +f dxbj,(x)(z — z')(x). So that 


J exedy26 70200406 998,0) + bj CODE x, Yyap(0) + dy (CD EC, 940) 


+a,(x)D¥(x, y)ag(y) = 0 
(66) 


where we use the shorthand notation D5';(z,x,y) := D5',(x,y) and 
similarly Dt¥(z;x,y) := Di (x,y)- 
Taking i€ {1,...,n}a € {1,...,p} and pe {1,... 


this as a matrix positivity condition 


,p+1}, we can write 


|z0 (67) 


/ day en ar | by) 


DP(x,y) Do(x,y) | La) 


where D2(x,y) is the (p+1)nx(p+1)n matrix-kernel with elements 
Dy i y), D?"(x, y) is the (p + 1)n x p matrix-kernel with rows labelled by 
piand columns labelled by 6 with elements DIB (x, »? and Dolx, y) isthe 
p x p decoherence matrix-kernel with elements pv (x,y). De ; describes 
the quantum back-reacting components of the drift. 

Equation (67) is equivalent to the condition that the 


((p + 1)n+ p) x ((p + 1)n+ p) matrix of operators 


ee a “6 


(68) 
Di Do 


be positive semi-definite. Here we are viewing the objects of 
Eq. (68) as matrix-kernels, so that for any position-dependent 
vector by (X), (Dz byi(x) = f dyDE' (x, y)b/ 0). 


6(Z,Z') 


Since we know D2 and Do must be positive semi-definite, we know 
from Schur decomposition that 


2D, > DP D;'pet (69) 


and 


(I — DyDo')D =0, (70) 
where D,' is the generalised inverse of Do. Furthermore, from Equa- 
tion (70), we see if Do vanishes, then clearly p™ must also vanish in 
order for (68) to be positive semi-definite. 


Observational trade-off in the presence of fields 

We can use the same methods to arrive at an observational trade-off 
using the field-theoretic version of the Cauchy-Schwartz inequality in 
Eq. (18). This arises from the positivity of 


Tr| / dzdz'dxdyA"”(z\z',x,y)O,(Z,2',X)p(z')O1 (z,z',y) >0 (71) 


for any local vector of CQ operators O,,(Z, 2’,x). We have to be careful, 
since (71) is not in general well defined since O,, may not be trace-class. 
We hence assume that we consider states p(z) and operators O,,(z, Z’, x) 
for which Eq. (71) is well defined. Since we are interested in getting an 
observational trade-off we expect this to always be the case for 
physical classical-quantum states p(z). 

We shall use Eq. (71) to arrive at a (pseudo) inner product on a 
vector of operators O, via 


(O,, On) = / dzdz'dxdyTr [Maiz'x, YOy(Xz)},(9)| (72) 


where ||O|| = \/ (O, O) = O due to Eq. (71). Technically this is not positive 
definite, but again, this will not worry us. Hence, so long as ||O>||#0, 
which holds due to the CQ inequality derived in the derivation of the 
Pawula theorem”, we again have a Cauchy-Schwartz inequality 

110,17 1105||? — (0, 02)|? 20. (73) 


Choosing O01, p(X) =6 eae) and 


On x= f dx’ Booz — 2'),00L,(x), one finds 


NOy?= | dedxdytr [Dgtz.x,yLeOde@LhO)] = (Do) (74a) 


\|O,|/? =2 i: dzdxdyTr [B(x )D4'(2:x, yOezZLI'Q)| (74) 


‘=(( il dxbi ‘cookin | 


(74c) 


(Oy, 0,)/? -| / dzdxTr [6° DF Z:x)Lg@alZ)L} 


Taking the limit b'(x) > 6(x,X)b'(x), we arrive at a local trade- 
off between diffusion, drift and total decoherence. In particular, 
using (74), the definitions of the expectation values of couplings 
defined in Eq. (33) and the fact that for back-reaction the 
expectation value of Do cannot vanish, we arrive at the observa- 
tional trade-off of Eq. (34) 


B(x) [2(Dp,(X,X)(Do) — (DP4)/?| Hx) 20 (75) 
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which we write in matrix form as 


2(D2(,X)) (Do) > (DE"(X)) (DP(X))". (76) 
It then follows directly from ved (76) that when the back-reaction 
is sourced by either D, ne or Dy components we can arrive at the 


observational trade-off i in terms ‘of the total drift 

a _ +t 

8(D2(X,X)) (Do) = (D} (X)) (Di (X))', (77) 

where in Eq. (77) recall that the definition of (Dt (x))" is given by Eq. (35) 
in the main body. 


A spatially averaged observational trade-off 

It is also useful to note that one can arrive at a spatially averaged 
observational trade-off which can be used to bound all of the elements 
of the diffusion matrix, not just its diagonals. Specifically, taking Eq. 
(74) with b'(x) = b' a constant, we arrive at the trade-off 


; 

8 / dvdy(D,0%,y))(Do)>( | dxb}x) )( | dxd{ ) F (78) 
where we define the expectation matrix 

Dxcey= | dzTr [D8 (Zs Db yOeZLt 0). (79) 


For the Newtonian limit discussed in the main body this bounds 
the diffusion in terms of the total mass of the particle 


2 


M 
80 
/ dxdy(D,(x, y)) = I6A (80) 


We can also arrive at a trade-off in terms of the eases New- 


tonian potential sourced by the masses by taking b‘(x) = =a In this 
case, we find the trade-off 
(Dx(x,9) DI(x) Di(x) \ 
3s / dxdy 2" _ (Do) > paw = fa <1 
|x — X||x — y| |x — x| |x — x| 


(81) 


which for the Newtonian limit gives a trade-off between the diffusion 
matrix and the effective Newtonian potential of the particle as sourced 
by its expectation value 


a 
i dxcdy Daron), [Fare _KOG=F (gay 
Ix—x||Ix—yl ~ 16d 16G*A ’ 


where we have defined the effective Newtonian potential 
(rue) 

as (®) = —Gf dx 
Newtonian limit of CQ theory 
In this section we motivate the Newtonian limit of gravity used in 
Section ‘Physical constraints on the classicality of gravity”. A fuller 
treatment can be found in*. We begin with classical General Relativity 
in the ADM formulation*®. To derive the Hamiltonian, we start from the 
3+1 split of the four metric 

= —(Nedty +g; (ax! +Nic dc) (dx! +N dc), (83) 
in which case, denoting @m,™%m aS canonical variables for the 
matter degrees of freedom, we can write the action for minimally 


coupled matter 


. Og; ; 
s= / dx(! od + 1, Be NH win), (84) 


where we are ignoring the boundary contributions to the action. Here, 


_[_ ¢ jipp, 16mG 1 a 

n=| tence P+ ce gi/2 (eugur’n 50 +H™, (85) 
3 

Hy = Bye +H, (86) 


are the Hamiltonian and momentum constraints and 77 is defined in 
terms of the extrinsic curvature tensor of constant ¢ surfaces, Kj, via 


(87) 


i= ero (Ki Key). 


It is useful to note that the matter densities H™, 1!" can be 
related to the matter stress-energy 7” via 


HM =,/EN’T, (88a) 


HE” = JeNT?. (88b) 
We here take the Newtonian limit of the gravitational field to be 
given by 


N= (1+ o) N=0, 8y= (1-7 


with ®(x) corresponding to the Newtonian potential. The choice of 
n= — © HO : is to ensure that 7 is canonically conjugate to M. A 
more detailed derivation starts from the full weak field metric, and 
gauge fixing of the shift can be found in*. Here, we find the effective 
action can be written 


(89) 


C2 ij 
a) 591 1! = - 16’, 


S= i d*x x(me oon Tm pes = Hew) (90) 
where the Newtonian Hamiltonian is given by 
yew =He t+ HO” +A (91) 
with 
5 (38 ens re) (92) 
8G 
the pure gravity Hamiltonian, and 
= / PxO(x)m(x) (93) 


is the interaction Hamiltonian, from which we see that non-relativistic 
matter couples to the Newtonian potential through its mass density 
m(x). In the case where we have the state of matter being described by a 
point particle 6(x—x(Q) of mass m the pure matter Hamiltonian would be 


=mc2 + (94) 


O' pip; 
(m) ae] 
Ho nm 
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The equations of motion for the gravitational degrees of freedom 
reads 


. 2 
= — Ame ” (95) 
2 
® 
ance (96 
©- gag — MOD : 


which, for 77@ = 0 yields the Newtonian solution for a stationary mass 
density. In a Louivile formulation the dynamics for the density 
PO, Te, x, pi) is given by 


Pp _ (m) pa op 3 _ 6p _ 
op 7 (Het Ho +P} 000 55+ | PAM 5’ (97) 


where the Hamiltonian and momentum constraints tell us that 
p(®, Me, x’, p;) should only have support over phase space degrees of 
freedom which satisfy the constraint 1 p(®, 77m) = 0. From Eq. (97) we 
can identify the classical drift associated to the back-reaction of the 
matter on the gravitational field from the MO) 585 term, so that 


Dp 


L ty X) = — MX). (98) 

In the classical-quantum case, we promote m(x) to an operator m. 
In this case Eq. (93) is the interaction Hamiltonian used in ref. 24 to 
study CQ gravity. We see from Eq. (97) that in any theory whose first 


moment reproduces the Newtonian back-reaction on average 


Tr [{H),0}] = / dx Tr ne an, a (99) 
must have a DP", given by 
(DPF 20) = — (r00), (100) 


from which the discussion at the beginning of the section “Physical 
constraints on the classicality of gravity” follows. 

Note that the present discussion is insensitive to the details of the 
theory provided it satisfies Poisson’s equation on average. None- 
theless, it’s interesting that when starting from the relativistic theories 
of refs. 11,39, we find that the weak field limit resembles the models of 
Tilloy and Di‘ osi®, which is discontinuous in ®, rather than the con- 
tinuous model of Di*osi”. This is because” allows for non-zero con- 
jugate momentum 7 with the kinetic energy of a different sign, while 
in ref. 81, the momentum is set to zero via the constraint equations and 
gauge fixing of the lapse. The discontinuity then arises because we are 
operating in the c > ~ limit, while in ref. 33, the discontinuity arises due 
to sourcing the Newtonian potential via a weak measurement process. 
We refer the reader to ref. 81 for details. 


Weak field CQ master equations 

Although the trade-off we derive does not depend on the parti- 
culars of the classical-quantum theory (provided it reproduces 
Newtonian gravity in the classical limit), we give two concrete 
examples for completeness. In ref. 30 we show that there are two 
classes of classical-quantum dynamics, one which is continuous in 
phase space, and one which has discrete jumps in phase space. 
We will give examples of each. Although they are the weak field 
limit of Oppenheim", it is worth stressing that taking the New- 
tonian limit entails certain coordinate choices and restrictions on 
the metric. For example, here, we have restricted ourselves to 
metrics of the form of Eq. (89). Any gauge fixing of General 
Relativity which is done before deriving the master equation, is 


generally not equivalent to taking the master equations of 
Oppenheim", and then taking the appropriate limit®. 


Continuous master equation 
For the class of master equations with continuous back-reaction, 
specifying that the first moment on average satisfies Eq. (100) is 
enough (up to drift terms which vanish under trace) to ensure the 
master Equation includes a term 


0e ad rR 6 60 . 
Bo = Hel), 0) — aH) f dPx| may gre + 2x) 


3 & 101 
+ [ &xd’y YS QnoTg@x) CX FoiX HO) ~~ 


+5 ;/ Pxd’Do(, x,x')([C0), [o, HV], 


where H, is the purely classical gravity Hamiltonian. We have 
taken the dynamics, i.e., the drift to be local in x, while we allow 
for the decoherence and diffusion terms to have some range. In 
this case, the evolution law is still local but correlations can be 
created’**. One can also add extra diffusion and decoherence into 
Eq. (101) which we do not consider here since it only leads to 
worse experimental bounds. However, adding additional diffu- 
sion in ® will generally be required in order to impose the con- 
straint m~0*. This master equation is close to the one 
considered in ref. 24, where the decoherence and diffusion ker- 
nels are chosen to be the ones discussed in the second example 
of Examples of Kernels saturating the decoherence diffusion 
coupling constants trade-off. This is the weak field limit of the 
simplest realisation in ref. 11. The case where the diffusion is 
spatially uncorrelated D,(x,y)=e(x—x’) a regulator which 
approaches a scalar delta function corresponds to the Newtonian 


limit of the diffusion term e(x— x’)(N(x)\/g0), {/g), o}}- 
Another natural diffusion kernel is 


D,(x,y)= — D,(1+ O(x’))A,, 6x, y), which can be understood as the 
Newtonian limit of the spatial diffeomorphism invariant kernel 
discussed in the section “Diffeomorphism invariant kernel”. 

One can supplement the kernel by some mechanism to control 
the diffusion. For example, a friction term such as 


F(Q)=Dp 5 f dxdx’dy{NX) /8@), (Vga ex — x’), HO)}e}. (102) 


In the weak field limit, this would add a term proportional to 


F(Q)* -/ dxdx’ (1o(X’)a) (103) 


foes (X) 


to the master equation of Eq. (101). Such a term would break Lorentz 
invariance since it sets a temperature scale, although this is not 
necessarily a deal breaker, since it is believed by many that quantum 
gravity is also likely to also have an anomaly. However, the friction 
term is a modification to D,(x), and if too large, could run afoul of 
precision tests of General Relativity, such as the orbital decay of binary 
pulsars. 


Discrete master equation 
An example of a discrete master equation satisfying Eq. (100) is 


20 = H_(®),0} — HHO”, e+ & - | dx 
le eed wee ee (1- 


(104) 


20) 
or) wena (x) 5 £09, a} |, 


with rt a dimensionless constant, and f(x) = yw. We have here inclu- 
ded hf and c to make it easier to compare with experiments. To 
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leading order, we could drop terms proportional to ®(x)/c’ in both the 
exponential and in N./g=1—2@/c? inside the integral over x. This 
gives 


de ; fond i yj 6 P 1 
ag 7 tHe(®), 0) — Hy”, o) + / dx [e" [Dare gix yews" Cx) — 5 (mlx), 0} + 


(105) 


These dynamical equations are supplemented with modified 
constraint equations as outlined in ref. 37. In any case, the trade-off in 
Eq. (31) is a statement independent of constraints and constraint pre- 
servation, at least in the weak field limit. 


Examples of Kernels saturating the decoherence diffusion cou- 
pling constants trade-off 
In this section, we give examples of kernels satisfying the decoherence 
diffusion coupling constant trade-off in Eq. (23). For any choice of 
kernel, we can compute the degree of diffusion it induces in precision 
mass measurements (see the section “Detecting gravitational diffu- 
sion”) and decoherence experiments (see the section “Decoherence 
rates”) which allows us to rule out certain kernels experimentally. 
Diffeomorphism invariance may single out having both Do(x,y) and 
D(x, y) approach delta functions at short distances’’, but other alter- 
natives may be possible. 

As a first example, we shall take the Lindbladian coupling to be 
Gaussian, taking 


9B 3 


“ A 
Do’, = — pr Ex(% y) (106) 


where g(x,y) is a normalised Gaussian distribution. The mass mo is a 
reference mass, and we shall take it equal to the mass of the nucleons 
which were considered in the section “Physical constraints on the 
classicality of gravity’, meanwhile A” is a coupling constant that 
determines the strength of the Lindbladian. 

It should be noted that with this choice of smearing function, 
the pure Lindbladian evolution appearing in Eq. (101) can be taken 
to resemble the Lindbladian part of spontaneous collapse 
models®”””"3"5, except here, there is no need to think about any ad- 
hoc field, nor think of the collapse as being a physical process. Rather, 
one necessarily gets decoherence of the wave function for free, via 
gravitationally induced decoherence!?>???>"6 

We now find the diffusion kernel D(x,y) using the coupling 
constants trade-off in (23). For simplicity, we shall assume the trade-off 
is saturated, and we will take the back-reaction to be local, so that 
(DP, (x,y) = (DP); (x) 6(x,y). In this case we find 


Dx, 9)= 5OPI 0 Ta ee OPO», (107) 


where gj/(x,y) is the kernel inverse of a normalised Gaussian 
distribution. 
It is shown in ref. 127, that the inverse distribution takes the form 


EV OY =F(X W(X), (108) 
where 
F(x,y)= Il 3 Ca(Fo)H on (= a , (109) 
i=1n=0 


n 2n 
and the limit N > ~ is taken. In Eq. (109) cq(r9) = 4 and d is the 
spatial dimension, so that x = (4%, X, ...,Xq). 


In total then, we arrive at the expression for the D2 which saturates 
the bound 


sory" on 2 Fal O F(x, WE v(x, YDP YN). (110) 


D(x, y) = 


If we further take the back-reaction that of the Newtonian limit in 
the section “Newtonian limit of CQ theory” 
(wey (x,y) = 4 16°" 67° 6(x, y) then we find the D, which saturates the 
bound is 


1 m2 
D,(x,y) = gaat (X, WEY). (111) 
(0) 


Another example of a Lindbladian coupling which is familiar in the 
literature is, 
ap 


(112) 
|x — =a 


DB (x x,y)= = 


For a single Lindblad operator, this is the coupling introduced in** 
used to reproduce a CQ master equation of gravity with a decoherence 
rate given by the Diosi-Penrose formula’* °°. Here we consider the 
special case where the x, y dependence of Do(x, y) is the same for all 
a, B which need not hold in general. The fact that it gives the same 
decoherence rate as Diosi-Penrose can be seen by plugging Eq. (112) 
into the classical-quantum master equation in Eq. (101). 

To invert the kernel in Eq. (112) we use the fact that 


ane 
~ an (pay) “OO 


from which one can immediately read of the generalised inverse 
(Do")ap(X,y) to be 


(113) 


0 es (114) 


(Do )ap(%¥) = V5(6(x,y)), 
where (Do ap are the matrix elements of the generalised inverse of Do. 
As a consequence, we find for this specific choice of kernel that the 
diffusion matrix saturating the coupling constants bound in Eq. (31) is 


ti 28a a : 
Di (X,Y) = 5 DIS) —_ V5(6X, D150), (115) 


where we have also assumed the back-reaction is local. Taking the 
back-reaction to further be that of the Newtonian limit of Eq. (101) 
(Dey (x) = 46°86" we find 


za 


D,(x,y)= ar W(X, y)). (116) 


This diffusion kernel is argued for on the grounds of having the 
fluctuations satisfy a Poisson equation, in ref. 24. 


Diffeomorphism invariant kernel 

Attempts to derive the constraint algebra of a generally covariant CQ 
theory” (and J. Oppenheim, “The constraints of a continuous realisa- 
tion of post-quantum-classical gravity", manuscript in preparation), 
motivates the spatial diffeomorphism invariant kernel 


Dix,x’)= — 5 DV gCONCog!a"'A,.6(x, x’), (a7) 
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where A, is the Laplace—-Beltrami operator. One can also consider the 
full 3 +1 kernel, via A” 6(x, x’)6(t, t’) along with the associated Green’s 
function of A“ but this is irrelevant for the Newtonian limit. It is 
however useful in removing the apparent asymmetry in the expres- 
sions below, since one must recall that the 6(x, x’) is a scalar in the first 
coordinate and a tensor density in the second, and likewise 6(t, t’) has 
an implicit lapse N(x’) in the second position. This kernel’s weak field 
limit is 
DIF (x, x)= — 3 DEVE" + Ox) Ay (x, X), (118) 
which is close to that of Eq. (116), but with a correction term that turns 
out to be important. 
Using Dy(x,x')= —4N./g6(x,x), the Lindbladian kernel in 
dimension d which saturates the trade-off for this diffusion kernel is 


1 
Do, ixt(%1X) = ip Ve ON) EGOS )GX, X), (119) 


with G(x, x’) the Green’s function for —A. It is a density in the x’ coor- 
dinate and a scalar in x. In the weak field limit, and to Oth order in ®(w), 
this gives the Diosi-Penrose kernel, Eq. (112). 

One could also consider the kernel 


DiM(x,x)= — 5 DV eweleMA.NC.x), (120) 
which in the weak field limit is 
DEE (x, x)= — 3 DEVE" A@(x)6(x, x’). (121) 


A comment on divergences 

The kernel examples given above give rise to divergent variance in the 
classical degrees of freedom, since in both cases the diffusion coeffi- 
cient diverges when evaluated at the same point D>(x,x). Though we 
do not have a general proof, this seems to be a general feature of the 
coupling constant trade-off: for the examples where we can compute 
the kernel inverse, at least one of D(x,x) and Do(x,x) diverge. A 
divergent D3(x,x) generally leads to a formally divergent classical 
energy production, whilst a divergent Lindbladian coupling Do(x, x) 
can lead to a divergent energy production in the matter degrees of 
freedom. The latter is related to the BPS problem®® of anomalous 
heating, although it isn’t necessarily equivalent since some kernels may 
diverge and be well-behaved from the point of view of energy pro- 
duction. This is not an issue from a conceptual point of view, since the 
only reason we expect energy to be conserved is due to Noether’s 
theorem, and Noether’s theorem doesn’t apply when the evolution 
isn’t unitary. 

In the standard BPS problem, energy production in open 
quantum field theory can be made small by renormalizing the 
Lindbladian coefficient Do(x, y) appearing in the master equation. 
Thus the problem is merely one akin to the hierarchy problem, 
where we are required to introduce another energy scale. How- 
ever, in the case of classical-quantum coupling, the coupling 
constant trade-off tells us that we cannot re-normalise Do(x, y) 
without affecting D (x,y). In particular, tuning the diagonals 
Do(x,x) to be arbitrarily small (large) has the effect of tuning 
D>2(x,x) to be arbitrarily (large) small: heuristically, one trades 
energy production in the classical system with energy production 
in the quantum system, and the relationship is fixed by the trade- 
off. On expectation, the total energy could be preserved, and the 
back-reaction can even slow down the flow of energy, but it’s 
unclear if this is enough. 


However, it is worth noting that while D,(x,x’) may appear 
to diverge at a single point as x — x’, when integrated over 
test functions, [ dxdx'D,(x,x)fEof(x’) is usually well-behaved. 
The kernels discussed above have this property. When it comes 
to physically relevant quantities, such as measuring the gravita- 
tional diffusion in tabletop experiments, it is the smeared well- 
behaved quantity that is physically relevant. However, in cos- 
mology, we typically take the constraint equation of General 
Relativity to be exactly satisfied at each point, and so one 
might imagine that 7%(x), and hence D2(x,x) is the relevant 
quantity (see the discussion in the section “Detecting gravita- 
tional diffusion”). However, one can set mmo=O via a gauge 
freedom*. In GR, the counterpart to the m2, kinetic term is 
Gyattm(x), and it’s perhaps worth noting that this quantity is not 
positive definite. It is also important to note that here, we have 
taken the weak field and c > ~ limit of General Relativity. At short 
distances when the diffusion becomes large, we expect this 
approximation to break down. One possible method of studying 
this problem rigorously would be through the regularisation 
properties of the classical-quantum path integral which we 
introduce in**”’. 


Decoherence rates 
In this section, we relate decoherence rates to Do, and also to the 
average (Do) = f dzTr[D*’(z;x, Lal XOLZ0)I- In particular, we shall 
show that the decoherence rate of a mass in superposition, is given by 
Eq. (131) in terms of the Lindblad operators and D® and canbe related 
to the quantity (Do). 

We consider the case of a quantum mass initially in a partially 
decohered superposition of state |L) and |R). We describe the 
quantum state using creation and annihilation operators 
w(x), W(x) on a Fock space, related to the usual momentum-based 
Fock operators as W(x)= [ dpe*a;. The mass density operator is 


defined via r(x) = my"(x)W(x), where m is the mass of the particle. 
We assume that the state remains well approximated by a state 
with fixed particle number, and the superposition can be taken to 
be distributions centred around x=x, and x= xp with total mass 
M, i.e. for a one-particle state we could’ take 


|L/R) = f Pxf 2 COP" (10). We will take them to be well sepa- 
rated so that f, (x)fp(x) = 0, and we take the separation distance to 
be larger than the scale of the non-locality in Do(x, y). Mathema- 


tically this means that (LIDS (2; x, YL p)La(X)|R) ~ 0 for any local 
operators Lq(x) and L(y). 


With this orthogonality condition, we can then (at least initially) 
consider the joint quantum classical state restricted to the two- 
dimensional Hilbert space of these two states so that the total 
quantum-classical system can be written as 

o(®, mo, t)= (ie, Ty,t) a(®,Tg,t) ), (122) 
a*(D,1g,t) Up(D, 1g,t) 
where u,(, 7, t) and up(®, 1, t) correspond to some subnormalised 
probability distribution over the classical states of the 
gravitational field. 

We define the total quantum state pg by integrating over the 

classical degrees of freedom 


Pa a / DODrpA®, To, t), (123) 


and we shall relate (Do) appearing in the trade-off to the 
decoherence rate of the off diagonals of pg. Integrating over 
the classical phase space in Eq. (9), one finds the following expression 
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for the evolution of po 


= = i DoDrg _ i[H(®, Ty), oA(®, Ny)] 
+ [ DoDI | dxdy[DPCO, 29:x, Ye, To, OL;0) 
i 
= ZDI(®, To: X Lg al), (®, Mo, Oh]. 
(124) 


In particular, one finds that the off-diagonals (L| ie |R) evolve in 
part according to the commutator, and in part due to the Lindbladian 
term 


J Poor | dxdy| L\DY(O, ro: La OO, ro, CLLO)IR) 


— 
2 


(125) 
DEO, to: X, YLMLZO La, C®, Te, HR) |. 


Care must be taken however, because both the quantum Hamil- 
tonian and the Lindbladian coupling constants depend on the classical 
degrees of freedom which are affected by the quantum degrees of 
freedom, and thus the evolution of the quantum system is non- 
Markovian in general. 

We shall now study the two terms appearing in Eq. (125) sepa- 
rately, starting with the first term. Since we assume that the state is 
well approximated by a state with fixed particle number then the 
contributions to the first term in Eq. (125) only come from 
terms where L,(x) and Lp(y) have the same number of creation and 
annihilation operators. To compute the expression, one 
commutes through the creation operators to act on the (L| bra, and 
picks up a term f,(x). Similarly, one commutes the annihilation 
operators to the act on the |R) ket, and picks up a term fp(y). As a 
consequence 


(LIDE®, 19; X, Lg XA, Ty, OL IR) ~ DE, to: X, Nf OF RW) = 0, 
(126) 


where the last equality follows from the fact that we are taking the 
masses to be well separated and the range of Do(x, y) is assumed to be 
much less than the separation between the masses. 

Hence, the evolution of the off-diagonals comes from the 
(off-diagonals) of the unitary evolution and the second term in Eq. 
(125), the so-called no-event term. The off-diagonals of the no-event 
term is 


1 ‘ 
-5 | Pov [ AxdyD$'(O, o;x, LILO), ©, Mo, )IR), 
(127) 


which is negative definite and acts to exponentially suppress the 
coherence. To see this, note that expanding out Q(®, 7, t) in terms of 
the approximate 2 dimensional Hilbert space 


A, Te, t) =U (®, Te, tL) (L| + Ug(®, To, O1R) (RI 


: (128) 
+a(Q, Tg, t)|L) (RI +a (Q, Tg, t)|R) (LI, 


and using the fact that the range of Do(x,y) is much less than the 
separation between the left and right masses, we can write the off- 
diagonals of the no-event term as 


1 A 
—5 | Poon yf ©, 19:x,9)((LILOMLa OIL) + (RILO)LQ(X)IR)) Lla(®, MR). 


(129) 


Equation (129) already expresses the fact that the off-diagonal 
terms will decay, and the particle will decohere at a rate determined by 
the integrand of Eq. (129). 

We can go slightly further when in the presence of a background 
Newtonian potential which is dominant, such as the Earth’s M,. The 
Earth’s background potential dominates over small fluctuations in © 
due to the particles and we can approximate Eq. (129) by 


- FO Ce. IN(LIL,ONL POV aCOIL) + (RILZO)LO)La(X)IR))(LIPgIR), 
(130) 


where the coupling DB (x, y) depends on the background Newtonian 
potential, but is otherwise phase-space independent. The result is to 
exponentially decrease the coherence (L|pq|R) with a rate A deter- 
mined by 


1 7 
= 5 | dxayogfx, y(LiLfOrLaCDIL) + (RLAONLLOIR)). 13H 


Let us now show that the (Dp) term appearing in the trade-off (34) 
is always less than (twice) this decoherence rate when in the presence 
of a background potential. Specifically, we show that 


(Do) = / DODiy / dxdyTr|DI(®, 11: x, YLZLq(X(®, To) <2A, 
(132) 


where we assume that we are in the presence of a background 
potential. To see this, we first expand out the CQ state in terms of 
Eq. (128) and use the fact that Dg has a range less than the separation of 
the masses. We then arrive at the following expression for the left-hand 
side of Eq. (132) 


[ Ponm. | dxdyDP®, 19:X, NLILZOLa IL) UL(, To, 0)+ (RIGO La OR UR®, To, 0), 
(133) 


In the presence of a background potential, this dominates the 
contribution to the decoherence and we are left with 


: dadyD%P(x, (LIL,OLa COIL) (LIPalL) + (RILZOVLa(0)IR)(RIpQIR)). 
(134) 


Due to the positivity of the CQ density matrix (L|p9|L) and (R\po|R) 
must both be positive. Furthermore, they must sum to one due to 
normalisation, from which Eq. (131) directly follows. 

It is also important to note that though J is the decoherence rate 
of a particle in superposition of L/R states, the bound (132) holds even 
for fully decohered masses in any mixture of |L)(L|, |R)(R| states. This 
can be seen directly from (133) which depends only on uw, Up. 


Decoherence rate examples 

In this section, we give an explicit example of a decoherence rate 
calculation. Importantly, we see that in general, the decoherence rate 
can depend on the probability density. This suggests that the terms 
appearing in the trade-off relation will need to depend on expectation 
values such as the expectation value of the mass at a point x, rather 
than a stronger bound in terms of the mass density. This is perhaps not 
surprising, since the decoherence rate itself can be thought of as an 
expectation value, being related to the average time it takes for off- 
diagonal elements to decay. In conclusion, this motivates us to 
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advocate for the volume of the wave packet to be included in the figure 
of merit in future interference experiments. 

We take the Newtonian limit master equation defined by Eq. (101). 
We ignore the unitary part of the evolution, since it will not directly 
contribute to the decoherence rate, and can be small for a free particle 
in superposition. From Eq. (101) we find the relevant evolution for the 
quantum state fo, obtained by integrating over the classical degrees of 
freedom to be 


1 ‘ i 
Fo=5 | Pxd’ yew. (biax),oq,MOMN). 135) 


We now compute the off-diagonal elements for a particle in 
superposition of orthogonal |L), |R) states 


0 
(1 Fela) = — [xd yDocx.»9 (em) — mpm, 0) ~ ma)LI2aIR 
(136) 


where m, (x) = (L|m(x)|L) and similarly for the right state. We see that 
the off-diagonals decay exponentially with a rate determined by 


A= / d’xd? yDo(x, ym, (x) — Ma(X)(M,(y) — Ma). (137) 


In the main body, and the previous subsection, we have assumed 
that the superposition of the particle is much less than the typical scale 
of Do(x, y). In this example, this means that we take the particles suf- 
ficiently separated so that we can approximate Do(x, y)m, (xX) mMp(y) = 0, 
in which case Eq. (137) is precisely the decoherence rate calculated in 
Eq. (131) with L(x) = m(x), as is to be expected. 

The natural decoherence kernel from the point of view of dif- 
feomorphism invariance is Do(x,y)=Do6(x,y), in which case for a 
particle of mass M, and uniform wave-packet volume V, Eq. (137) gives 
A=2DoM/V. 

As another example, we can take Do(x, x’) to be the Diosi-Penrose 
decoherence kernel defined via Do(x,y)= ep so that the off- 
diagonals decay exponentially with a rate proportional to the 
Diosi-Penrose decoherence rate 


A= / dxd?y (138) 


D 
Ben ye) — Mag) (IM) — MARL). 


-Jl 


In this example, taking the superposition to be sufficiently sepa- 
rated means that we are approximating are =O in comparison with 


the rest of the terms appearing in Eq. (138). We are then left with 


A= / Pxd?y 


(139) 


D 
ie py UML) + Ma(X)IMa(Y)), 


Jl 


which for spherical distributions of radius R and total mass M is pro- 
portional to the average gravitational self-energy of each mass 
distribution A= 0M", 

For a composite particle of mass M, made up of N constituents 
each of radius R, the mass density will be represented by a sum over all 


of the particles m(x) = mx). The decoherence rate is given by 


D 
A= [ @xdy oo $s my, (20), (9) + S> Mp, (CMe, 1) : 
ij ij 
(140) 
Since the cross terms involving i,j are suppressed by a factor of 


inter-atomic scales, to leading order the contribution to the deco- 
herence rate is lower bounded by the i= component of the sums in Eq. 


(140), which gives an extra factor of N relative to the single particle 
case A= S20NM" 

Both Eqs. (138) and (139) depend on the probability density of the 
mass. In particular, taking the probability density to be arbitrarily 
peaked, one finds that the decoherence rate also diverges. This has to 
be the case: recall from the section “A Trade-off between decoherence 
and diffusion” that if one considers a particle in a superposition of two 
arbitrarily peaked probability densities, then there can be an arbitrarily 
large response in the Newtonian potential around those points. As a 
consequence, for such states, the decoherence must occur arbitrarily 
fast, or there must be an arbitrarily large amount of diffusion to cover 
up the back-reaction and maintain coherence. For the continuous 
master equation, such as that of Eq. (101) this diffusion must also occur 
throughout space, although it can depend on the gravitational degrees 
of freedom. Since divergent energy production throughout space is 
clearly unphysical, it must be the case that the decoherence rate must 
also depend on the expected mass density, as is the case for this 
example. This argument allows us to rule out continuous master 
equations that have pure Lindbladian terms that predict decoherence 
rates which remain finite as the mass density becomes arbitrarily 
peaked since the coupling constant trade-off will demand that an 
infinite amount of diffusion is required to cover up the back-reaction 
and maintain coherence. This is the case for the class of models with 
CSL-type Lindbladian couplings which are phase space independent, 
for example in Eq. (106). 


Detecting gravitational diffusion 

In this section, we show how the diffusion induced on the Newtonian 
potential can be measured experimentally. metric leads to observable 
effects, such as variations in the accelerations involved in torsion 
experiments and stochastic wave production in cosmology. As shown 
in the main body of the text, in the non-relativistic limit, c > ©, the CQ 
dynamics can be approximated by sourcing the Newtonian potential 
by arandom mass term, and in order to maintain the coherence of any 
mass is superposition, there must be noise in the Newtonian potential 
such that we cannot tell which element of the superposition the par- 
ticle will be in 


V0 =4nG[m(x, t)+u(®, MJ(x, o)], (141) 


with 


mlI(X, 0] = 0, EmlW(x, QUI, C')| = 2(Dy(x,y, ))6(t, 0), (142) 


where (D,(x,y,®)) := Tr[D5"(x,y, Dp )L,(X)PL)(y)] and p is the quan- 
tum state for the decohered mass density. The diffusion coefficient in 
Eq. (142) is chosen in order for the dynamics to have the same 
moments as the CQ master equation (4). The solution to Eq. (141), 
having absorbed u into / is given by 


(t,x) eee c | d3x’ [mx’, t) ae Ne i) , (143) 


where the statistics of J are described by Eq. (142). A formal treatment 
of solutions to non-linear stochastic integrals of the form Eq. (141) can 
be found in ref. 82. 

One can also verify this behaviour in specific cases. In the con- 
tinuous model of the section “Continuous master equation”, the noise 
is taken to be Gaussian, and this, as well as the evolution of the 
quantum state, is what determines the diffusion in Eq. (143). For the 
class of discrete models, the higher order moments such as 
EnV, O/(y, UY(z, t’)] are suppressed by an order parameter””°”” and 
whenever this is true we expect we can approximate the dynamics of 
the Newtonian potential by a Gaussian process. Whether this is the 
case or not, it is the second-order moment that enters into our 
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discussion of the variance here. As such, for minimally coupled the- 
ories, the Newtonian potential will appear to be sourced by a random 
mass distribution. 

In the discrete case, a precise understanding of the effects of the 
diffusion beyond the Gaussian approximation involves solving the full 
classical-quantum dynamics, perhaps using the methods of Oppen- 
heim et al.°. In Eq. (141) we are also taking the time scale of the dif- 
fusion to be faster than the dynamics of the matter distribution. 
Likewise for the decoherence—we showed in the “Methods” subsection 
“Decoherence rates” for continuous models the evolution of the 
quantum state acts to decohere it into a mass density eigenbasis m(x). 
One could of course also include the quantum state evolution in a 
simulation of full CQ dynamics, but this is beyond the scope of the 
current work. 

Note that if the moments of the noise process are Galilean 
invariant, then the theory given by Eq. (141) is Galilean invariant as we 
would expect in the Newtonian limit. One can further derive Eq. (141) 
from manifestly diffeomorphism invariant theories”, such as those of 
refs. 39,131. 


Table-top experiments 
In this section we estimate the variation in force which would be seen 
in table-top experiments which bounds the diffusion of classical the- 
ories of gravity from above, giving a squeezed bound on D, due to 
lower bounds on diffusion arising from coherence experiments. We do 
this for dynamics in Eq. (141), but the methodology is general and 
could also be used in a full simulation of CQ dynamics. 

The variation in force induced on a composite mass is found via 


> 


P= -| Pxm(x)VO. (144) 


Using the solution in Eq. (143), the total force can be written 


B= 2) 


se Ime) JO, 0) 


F ot= -G / d?xd?x' mex) (145) 

In reality, we measure time-averaged force by measuring time- 
averaged accelerations over the time resolution of the experiment 
ATS il dtF,,.,. The total variation in the force’s time-averaged mag- 
nitude. The full covariance matrix for various kernels in the Newtonian 
limit is ues in J. . Oppenheim and A. Russo, manuscript in prepara- 
tion), o7 := F tot * oF tot can be written as 


1 ; ¥-X)G-y - 
oz= ape / xd? yd?x'd? ym) Seer eae (Dx(x', y', ®)). 


(146) 


We shall use Eq. (146) to provide an upper bound on coupling 
constants of CQ theories for different choices of kernels D(x’, y’, ®). 
Given a choice of functional form of the kernel, all that remains is the 
strength of the diffusion coupling, which for the translation invariant 
kernels we consider here takes the form of a single coupling constant 
D>. We take D> to be a dimension-full quantity with units kg’ sm? which 
characterises the rate of diffusion for the conjugate momenta of the 
Newtonian potential. 

For a composite mass, we can approximate the mass density by 
summing over N individual atoms of mass density m((x), m(x) = ymx). 
The total force is the Ben by F tot = >i F ;, where F ; is the force on 
each individual atom F ;= — : yazmeoN Oe, and the total variation 
of force is then 02 = pF — EDF’ 

In general, the squeeze will depend on the functional choice of 
D2(x, y, ®) on the Newtonian potential. As mentioned in the main body, 
in the presence of a large background potential ®,, such as that of the 
Earth’s, we will often be able to approximate D(x, y, D) = D2(x, y, Dp). 


This is true for the kernels with functional dependence of the form 
Dz ~ 0", D,~V ®, though the approximation does not hold for all ker- 
nels, for example, D, ~ V’® which creates diffusion only where there is 
the mass density. We hereby shall only consider diffusion kernels 
D>(x, y, Pp) where the background potential is dominant, leaving more 
general considerations for future work. 

For local translation invariant dynamics for which the background 
Newtonian potential is dominant, for example, D,~®", we have 
(D2(x, y, Dp)) = (D2(Mp))6(x, y) and we arrive at the expression for the 
total variation in time-averaged force 


-*)-G- 
Ix — x’ Ply — x’ 


ul d?xd? ‘yd? x’ mj(x)m; on Ving ®,)). 


> 
(147) 


To leading order, the integral in Eq. (147) is dominated by the self 
variation term where i=/, since nuclear scales 10°° m dominate over 
inter-atomic scales 10 m, so that ODay iFl~ Di [F?]. Approximat- 
ing the mass density of the atoms as coming from their nucleus, and 
taking them to be spheres of constant density p with radius ry and mass 
my, We find that the integral in Eq. (147) is approximately 


NG? 
on Ne eN | ax (Dy(,)). 


(148) 


For the class of continuous dynamics (D2(®,)) = D2(®,), since the 
diffusion is not associated with any Lindblad operators. If there is noise 
everywhere throughout space, then the integral in Eq. (148) diverges 
and gives evidence that continuous CQ theories with noise everywhere 
should be ruled out. 

As such, we expect that continuous CQ theory must contain non- 
linear terms proportional to the Newtonian potential appearing in Eq. 
(141), in which case we can approximate { dx’D, by V,D2 where Vj is 
the volume of the region over which the background Newtonian 
potential is significant. In total then, we find for continuous local CQ 
dynamics 


2 42 -2 
alex D,NG' p*ryVy, (149) 
: AT 
From this, we can calculate D> in terms of the total variance of the 
acceleration 03 = = mz, to geta lower bound 


< NAT (150) 
VG 

Standard Cavendish type classical torsion experiments measure 
accelerations of the order 10’ ms”, and we can take the time over 
which the acceleration is averaged to be that of minutes AT -~10*s,soa 
very conservative bound is o,~107’ ms”, whilst N will be N~ 107° and 
ry~1075 m. We take the background Newtonian potential to be that of 
the Earth and we (conservatively) take V, to be V, ~r2h~10” m? 
where rg is the Earth’s radius and h is the atmospheric height. We see 
that this bounds D, from above by D, < 10 kg? sm*°. 

On the other hand, D, is bounded from below from interferometry 
experiments which bound the decoherence rate. From Eq. (137) and the 
coupling constant trade-off, for the kernel D (x, y) =D26(x, y) we see 
(ignoring constant factors) that the decoherence rate is found to be 


NMi 


151 
V,D,’ 6) 


where M, is the mass of a composite particle in the interferometry 
experiment, which is made up of N, particles, each with volume V,. This 
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gives rise to the squeeze 


ORNryAT | 
VC 


, NaMi 
25°V A 


(152) 


Using the numbers from ref. 59, with M,~10* kg, N,~ 103, and 
V,~ 10°10 107 m? =10~” m?, A-10's™ we find that D2 >10° kg? sm™. 
This suggests that the D(x, y) = D26(x, y) kernel for classical gravity is 
already ruled out by experiment. 

For the local discrete models, such as that of Eq. (104), the theory 
is less constrained due to the dependence of the diffusion on the mass 
density. In this case (D,(®,)) = #Dy(,)m(), where the factors of 
Planck length and Planck mass are to ensure that D2(®,) has the 
required units. We arrive at the upper bound for D 


o2Nr} ATmp 


2D. 
myG@'lp . 


(153) 


Meanwhile, from Eq. (131), and coupling constant trade-off (30) 


the decoherence rate for local discreet jumping models goes as 

Aw sare, which gives rise to the lower bound for D>. From this, we 
pY2 

arrive at the squeeze 


OANrWAT | [pD, | My 


>P2s i, (154 
myG Mp A 


and plugging in the numbers we find the bound given by Eq. (46) which 
gives rise to the squeeze for local discrete mod- 
els 10-'kgs > “.D, >10~* kgs. 

We can also consider other diffusion kernels, for example, that of 
Eq. (118). In this case, for continuous dynamics, we have that 
(D(x, y)) = — ED,(,)V6(x, y). The Lindbladian kernel saturating the 
coupling constants trade-off at zeroeth order in (x), is the 
Diosi-Penrose kernel Do(x,y, ®,) = Pols), as we saw in the section 
“Examples of Kernels saturating the decoherence diffusion coupling 
constants trade-off’. Approximating the masses as spheres of constant 
density we find from a substitution of the kernel into Eq. (146) that the 
variation in time-averaged force is given by 


a) 
o2~ Ns (155) 
ATry 
We therefore find a lower bound for D> in terms of the variation in 
acceleration 


2 
_ ATIRORNF} 


Dys 3 (156) 


, 


which for classical torsion experiments o, ~ 107’ ms”, T~ 107s, N~ 107° 
and ry~10 m gives Dot <10~°kg m-!. On the other hand, for this 
kernel the decoherence rate can be calculated via Eq. (139) 


2 
soe. (157) 
(DR, 
which gives the squeeze on D 
ATo2Nr3 N,M 

oA Shp R oa (158) 
For the numbers used in the main body of the 
text?°M,~ 10 kg, Ny~ 107, R,- V2 =10"m,A-10's, this _ yields 


Dt >10-* kgm and so this model is not ruled out by experiment. 


In general then, we expect that by simulating full CQ dynamics 
satisfying the decoherence diffusion trade-off we will be able to 
squeeze Dz from above and below. We bound D, from above by 
studying the effects of diffusion on gravitational experiments, and we 
bound D, from below using the coupling constant trade-off and 
coherence experiments lower bounding the decoherence rate. As we 
have seen in this section, it appears that classes of continuous CQ 
hybrid theories of gravity, including models without spatial correla- 
tions, are already experimentally ruled out, whilst others, such as the 
kernels in subsection “Diffeomorphism invariant kernel” require 
stronger bounds from both gravitational and coherence experiments. 
We have been very conservative in our estimates, and so we expect a 
more thourough analysis will tighten the bounds by orders of 
magnitute. 
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